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FOREWORD 


The investigation described herein, which was conducted by the Liquid 
Rocket Operations, Propulsion Division, of Aerojet-General Corporation, 
Sacramento, California, was performed under Contract NAS 3-7977. The work 
was done under the management of the NASA Project Manager, Mr. W. A. Rostaf inski 
of the Chemical Rocket Division, NASA Lewis Research Center. 
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ABSTRACT 


This six-task effort was initiated in March 1967 and completed in 
March 1969. Its purpose was to formulate and experimentally verify an 
analytical computer model for the transient performance of a two-speed 
inducer system, hydraulically-designed for operation in liquid hydrogen. 

A two-speed inducer system was designed, fabricated, assembled, and tested. 

A hydraulic turbine, powered by the flow from a rotor connected to a high- 
speed drive, drives the low-speed inducer. The high-speed rotor blading is 
mounted between the low-speed inducer and the hydraulic turbine blading with 
essentially all of the delivered flow passing through each row of blades. 

The system was designed to operate at a flow rate of 4900 gpm with a 25 ft 
NPSH and an over-all cavitating head rise of 5453 ft; the high-speed rotor to 
operate at 44,500 rpm and the low-speed inducer at 18,800 rpm. This system 
could potentially reduce the required NPSH of a liquid hydrogen engine. The 
test results indicate that the system is dynamically stable, but additional 
research is needed to improve the cavitation performance of the high-speed 
rotor. The computer model was successfully developed. 


xi 



I. 


SUMMARY 


The objective of the six-task Inducer Dynamics - Full-Flow, Full- 
Admission Hydraulic Turbine Drive Program was to formulate an analytical com- 
puter model for the transient response of a two-speed inducer system which 
was hydraulically designed for operation in liquid hydrogen. This model was 
verified experimentally, using water as the test fluid. 

Summaries of each task follow along with the major conclusions and 
recommendations . 

A. TASK I — HYDRODYNAMIC DESIGN AND MECHANICAL LAYOUT 

A two-speed inducer system was designed with the low- speed inducer 
being driven by a hydraulic turbine. The power source for this turbine is 
the flow supplied by a rotor connected to a high-speed drive. The high-speed 
rotor blading is mounted between the low-speed inducer and the turbine blading 
with essentially all of the delivered flow passing through each row of blades. 
The design parameters are as follows: 

NPSH - 25 ft (LH 2 ) 

Flow Rate - 4900 gpm 

High-Speed Shaft Rotational Speed - 44,500 rpm 
Low-Speed Shaft Rotational Speed - 18,800 rpm 
Head Rise (Cavitating) - 5435 ft 
Efficiency (Cavitating) - 61.3% 

A design layout for an 8500 rpm high-speed shaft water test dem- 
onstration unit was completed. The low-speed shaft of this assembly is sup- 
ported by hydrostatic bearings and can be braked by means of a Prony Brake 
to alter the normal operating speed ratio. This Prony Brake also is designed 
to be utilized as a dynamometric device for measuring torque. In addition, 
design layouts were completed of the assemblies for driving the high-speed 
shafts with an electric motor for steady-state testing and with a gaseous 
nitrogen turbine for transient testing. The test unit was designed to permit 
its conversion to liquid hydrogen operation without any need for major 
modif icat ions . 

B. TASK II— FORMULATION AND DEMONSTRATION OF COMPUTER SIMULATION 

A model was formulated which describes the design and off-design 
pumping system performance. Along with this model, an associated digital 
computer FORTRAN IV program was formulated for use with the IBM 360 Model 65 
computer. This program also is compatible with the Model 67 computer. 

The steady-state and transient dynamics of the flow of fluids in 
lines that are external to the pumping system are described using the method 
of characteristics to evaluate an elastic column "waterhammer 11 model. An 
adequate number of test cases were run to verify the capability of the model 
and the computer program. 



C. TASK HI — PARAMETRIC STUDIES 

The effects from changes in the following were evaluated in 
parametric studies: 

- High-speed rotor speed ramp rate 

- Mass moment of inertia (low-speed turbine 

and hydraulic turbine) 

- Pump suction pressure 

- Suction line geometry 

- Discharge line geometry 

- Discharge valve position 

- Fluid acoustic velocity 

- Pumping system size 

D. TASK IV--DETAILED DESIGN 

The detailed design effort was minimized by using standard MS, AN, 
NAS, or commercial parts wherever practical. An existing NERVA power trans- 
mission was modified slightly and used in the design. 

E. TASK V — FABRICATION AND ASSEMBLY 

Three blade rows were fabricated and only minor modifications are 
needed to test this hardware in liquid hydrogen instead of water. Fabrication 
and first assembly costs were kept relatively low by utilizing a single vendor 
for all procurement, complete fabrication, and initial assembly. As a result, 
only minimal vendor coordination was needed and minor assembly problems were 
readily foreseen and corrected. 

F. TASK VI— TEST AND COMPUTER MODEL REFINEMENT 

The data needed as input to the computer model was obtained from 
the testing conducted with water. The test unit operated satisfactorily at 
flows from shut-off to 125% of the design value. The low-speed shaft remained 
hydraulically- locked with the high-speed shaft over the full range of flow and 
the Prony Brake made it possible to obtain data at the extreme off-design 
speed ratios. The large, shroud-mounted, hydrostatic thrust and journal bear- 
ings were successfully operated during the complete test program. Low-speed 
inducer head-flow performance was below the design value, probably as a result 
of choking in the blade row at the entrance to the airfoil section. The high- 
speed rotor cavitation performance was lower than expected, but the cause is 
unknown because there is no design criteria for high stagger angle blading. 

The following is a comparison between the over-all pumping system 
water design point with the over-all pumping system water test results: 
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Design 

Test 

Flow Rate 

935 gpm 

935 gpm 

High-Shaft Speed 

8500 rpm 

8500 rpm 

Low- Shaft Speed 

3590 rpm 

3100 rpm 

Head Rise 

200 ft 

150 f t 

Efficiency 

61.3% 

43% 

The transient tests were successful 

in that the speed ramps 

were easi] 


trolled and the low-speed shaft followed the high-speed shaft except during 
extreme cavitation of the high-speed rotor. The computer model prediction 
agreed with the start transient test data for the critical parameters of 
speed, flow, and discharge pressure. Also, it was found that predicted tran- 
sient was sensitive to the hydraulic turbine torque characteristics. 

G. MAJOR CONCLUSIONS AND RECOMMENDATIONS 


Based upon the results of this program, the following are the 
major conclusions and recommendations associated with the two-speed boost 
pump (viz., full-flow, hydraulic, turbine-driven inducer) investigated. 

1. Appropriate evaluation must be made of the low-speed inducer 
and hydraulic turbine torque/speed matching characteristics to ensure that 
the pumping system is stable and hydraulically-locked . 


2. The more conventional high-speed rotor blade geometries of 
the lower head rise (1500 ft to 2500 ft) designs are more easily mated to the 
downstream main-stage pump. 

3. A tandem blade row low-speed inducer with a high inlet flow 
coefficient should be considered for future designs. Then, the low-speed inducer, 
which is a necessarily high head coefficient machine would exhibit increased 
performance predictability as well as increased inducer efficiency. 

4. The mean radius from the low-speed inducer inlet to the 
hydraulic turbine inlet should be maintained constant if at all possible 
although some increase in inducer radius is permissible. Whatever change in the 
mean radius is necessary to mate the flow channel with the downstream main- 
stage pump should be made through the hydraulic turbine (radial end flow type 

or mixed flow type). 

5. Design criteria should be developed to establish typical 
high-speed rotor blading cavitation performance. This blading should have 
minimum fluid angles of 80-degrees at the tip and fluid loading producing a 
diffusion factor of 0.50. 

6. This type of boost pump should be used in rocket engine pump- 
ing systems that permit the main pump to be mounted close to the tank. The 
main-stage pump will operate two to three times faster than the inducer result- 
ing in a compact pump with minimal over-all weight. 
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7. If the off-design characteristics of each blade row are known 
or can be predicted, the computer model can be used successfully to study the 
start characteristics of a particular engine system. 

II. INTRODUCTION 


The work performed by the Liquid Rocket Operations of the Aerojet-General 
Corporation, Sacramento, Colifornia, to successfully accomplish the objectives 
of the Inducer Dynamics - Full-Flow, Full-Admission Hydraulic Turbine Drive 
Program are delineated herein, This effort was conducted under Contract 
NAS 3-7977 and the cognizance of the NASA Lewis Research Center. 


The following six-task contractual effort was initiated during March 1967 
and completed in March 1969. 


Task I Hydrodynamic Design and Mechanical Layout of Hydraulic- 
Turbine-Driven Low-Speed Inducer/High-Speed Inducer pump 
System 

Task II Formulation and Demonstration of Computer Simulation to 
Predict Transient and Steady-State Behavior 

Task III Parametric Studies 


Task IV Detailed Design 

Task V Fabrication 

Task VI Pump System Tests to Determine Performance and to Develop 
the Computer Simulation 


The full-flow, hydraulic- turbine-driven inducer concept, which is shown 
on Figure No. 1, incorporates a two-speed power transmission and pumping 
system. The low-speed inducer is driven by a hydraulic turbine, the power 
source for which is the flow supplied by a rotor connected to a high-speed 
drive. The high-speed rotor blading is mounted between the low-speed inducer 
and the turbine blading with essentially all of the delivered flow passing 
through each row of blades. In a rocket engine feed system, the high-speed 
rotor would be attached directly to the main pump rotor. The operation of the 
low-speed inducer at speeds below those of the main pump results in feed system 
NPSH requirements that are less than those of the directly-coupled inducer 
systems. 


The full-flow, hydraulic- turbine-driven inducer concept has been con- 
sidered for the advanced NERVA nuclear rocket application (Ref. 1 and 2). 
Designs of this type also are being considered for the Aerospike Nozzle Concept 
in the Advanced Cryogenic Rocket Engine Program (Ref, 3). 


Aerojet-General has previously designed, built, and tested hydraulic- 
turbine-driven inducers wherein the hydraulic turbine operated using only a 
portion of the delivered flow (Ref, 4), which is supplied by means of 


4 




Figure 1. - Full-Flow, Hydraulic- Turbine-Driven- Inducer (Axial Section), 








recirculation from the main stage pump. The test program associated with 
these pumps consisted of steady-state tests at Aerojet-General and transient 
tests conducted by the Air Force Rocket Propulsion Laboratory (Ref. 5) using 
the same hardware. 

Testing of the inducer during Task VI was accomplished using water as 
the test fluid; however, the unit is designed to comply with the contract 
specification in that it can be tested in liquid hydrogen without extensive 
modifications. These Task VI pump system tests involved both steady-state 
and transient testing. The steady-state operating range of the pump was from 
80% to 120% of the nominal design flow coefficient (flow rate/shaft speed) 
with the high-speed shaft operating at speeds of up to 8500 rpm. 


The contractual requirements for 

are: 

Fluid : 

Flow Rate: 

Minimum Head Rise: 

Net Positive Suction Head: 
Off-Design Performance: 


the test system hydraulic design point 


Liquid Hydrogen 

4900 gpm to 11,400 gpm 

5000 ft (at hydraulic turbine exit) 

25 ft 

To 120% of design flow coefficient 
(flow rate/high-speed shaft speed) 
without cavitation in the high-speed 
rotor when operating at design speed 


The ensuing sections of this report detail the technical effort for 
each task as well as present the conclusions and recommendations resulting 
from this effort. 


6 



III. TASK I - HYDRODYNAMIC DESIGN AND MECHANICAL LAYOUT 


The design of the demonstration test unit is detailed in this section. 
This design effort was divided into the following three subtasks: 

- Hydraulic design of a hydrogen test unit. 

Design of a water test unit driven by an electric 
motor to be used for the steady-state performance 
testing. 

Design of a water test unit driven by a gas turbine to be 
used for the transient performance testing. 


A modified MARK III, MOD 4 NERVA Experimental Engine turbopump turbine 
was used to drive the water test unit for the transient testing. Only the 
first stage of this two-stage turbine was used with gaseous nitrogen at 
ambient temperature as the driving fluid. In addition, the turbopump/bearing 
housing was modified to accept grease-packed bearings. 


Adequate design work was accomplished to ensure that the test unit 
could be converted to liquid hydrogen testing without need for any extensive 
modifications. The single-stage turbine would be utilized in combination with 
the unmodified, liquid-hydrogen-lubricated bearing system. However, the 
maximum operating speed of the turbine is approximately two-thirds that of 
the selected hydrogen hydraulic design speed as a result of turbine stress 
and vibration limitations. 

A. DESIGN POINT SPECIFICATION 
1. Over-All Performance 


The following values were selected as a design point: 

Flow rate; 4900 gpm 

High-Speed Shaft Speed: 44,500 rpm 

Head Rise (Cavitating) : 5435 ft 

The contractually-required minimum flow rate was selected to 
minimize liquid hydrogen consumption for potential hydrogen test applications. 
In addition, the selection of this flow rate permits partial-speed testing of 
the full-scale hardware in available test facilities. 

An analysis of existing as well as planned hydrogen pumps 
was made to select the high-speed shaft rotational speed. A preliminary 
value of 30,000 suction specific speed was selected for the low-speed inducer 
based upon the use of working fluids without thermodynamic vapor head 
depression characteristics. However, the contractually-specified NPSH of 
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25 ft is for hydrogen, which exhibits head depression characteristics; 
therefore, a correction of 80 ft was applied to obtain a value of 105 ft. 
This latter value was used for the calculations of suction specific speed. 
Both the 30,000 suction specific speed and 80 ft liquid hydrogen head 
depression values are representative of current liquid hydrogen inducer 
design values. The 80 ft value was used in designing the M-l engine fuel 
pump (Ref. 6). 


An equivalent low-speed inducer speed was calculated for each 
of the pumps shown on Table I using the above indicated values and flow rates. 

The ratio of turbopump speed to low-speed inducer speed also was calculated 
using the actual operating speeds of the turbopumps* Then, the various pumps 
were categorized by speed ratio as shown on Table II. A speed ratio of 3.15 
(corresponding to that obtained in a large NERVA turbopump study) was selected 
as being representative of current turbopump applications. The high-speed 
shaft design rotational speed was calculated using the selected design flow 
rate and was based upon the 3.15 speed ratio, the low-speed inducer suction 
specific speed of 30,000, and the 105 ft NPSH values. 

Further analysis, subsequently discussed, indicates that 
the full-flow inducer device would not operate in a stable condition (low- 
speed inducer speed hydraulically-locked to the high-speed rotor speed) at 
the originally selected design speed ratio of 3.15. Reduced design speed 
ratios were required to obtain stable operation. Rather than alter the 
selected high-speed shaft speed, the low-speed inducer speed was increased. 

This corresponds to a suction specific speed increase from 30,000 to 40,000, 
and a speed ratio decrease of from 3.15 to 2.37. Analysis indicates that 
stable operation is obtained at this condition. 

An evaluation of typical boost pumping system head rise 
requirements was made to ascertain a suitable design head rise value. 

Table III summarizes the results. The contractual minimum requirement of 
5000 ft is a reasonable intermediate value for the various possible applications. 
A non-cavitating head rise value of 5600 ft was selected to establish blade 
mean line design parameters. The additional 600 ft of head was added to allow 
for low-speed inducer head loss resulting from cavitation; the differences 
between the preliminary mean line and the final blade element, integrated 
head rise; and prediction uncertainties. The final cavitating head rise 
obtained from blade element predictions is 5435 ft (see Section III,B) . 

2. Stage Work Split Selection 

The following factors were considered in selecting the stage 

work split: 


- stage efficiency and blade loading 

- interstage power/ speed matching characteristics 

- stage cavitation 
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SHAFT SPEED RATIO FOR DUAL PUMP CONCEPT APPLIED TO CURRENT 
AND FUTURE LIQUID HYDROGEN ROCKET ENGINES 
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(2) Ratio main pump rotor speed to low speed inducer rotational speed. 



CATEGORIES OF TURBOPUMPS ADAPTABLE TO DUAL-SPEED PUMP CONCEPTS 
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TYPICAL MAIN STAGE HIGH-SPEED PUMP PARAMETERS 
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Analyses were conducted using blade mean line calculations. 
Design was restricted to staging, which consisted of rotors only because pre- 
liminary analyses indicated that stators offered no significant design point 
performance advantage while they reduced off-design cavitation performance 
capability. The presence of stators in the staging also complicates the mech- 
anical design (see Section III,C). To further simplify design and fabrica- 
tion of the test unit, the high-speed rotor blade and the hydraulic turbine 
blade annulus dimensions were made identical (i.e., blade height and mean line 
radius were equal and constant). Preliminary analyses indicated that reasonable 
designs could be achieved with such geometries. 

a. Stage Efficiency and Blade Loading 


For the low-speed inducer stage, available test data 
for inducers were correlated by plotting loss coefficient and relative 
velocity retardation factor versus theoretical head coefficient as shown on 
Figure No. 2. The parameters are defined as: 


co 


(H. 


ideal 




/2g) 


co = loss coefficient 


^ideal 


ideal (isentropic) head rise, ft 


H 


actual 


actual head rise, ft. 



g - 

Ul _ 

ideal 


inlet mean (rms radius) relative velocity, 
ft/sec 

2 

acceleration due to gravity, ft/sec 

H.j . /U 0 , - Cu 0 /U„ - ideal head coefficient 

ideal 2m/ g 2 2m 


Cu 


2m 


U 2m 

V W lm 


w 


2m 


exit mean (rms radius) tangential absolute 
velocity, ft /sec 

exit mean rotor velocity, ft/sec 
relative velocity retardation factor 
exit mean relative velocity, ft/sec 


- inlet mean relative velocity, ft/sec 

The relative velocity factor was utilized, rather than 
the diffusion factor, D, (Ref. 7) to simplify calculations. Because of the 
low coefficients, high blade stagger angles, and larger secondary flows in 
inducers, there appeared to be no advantage in using the more rigorous two- 
dimensional correlation. Data for helical inducers (Ref, 8, 9, 10 and 11) 
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and inducers with larger head rise, including single blade row (Ref. 12) and 
tandem blade row (Ref. 13 and 14) designs are shown. The cross-hatched line 
represents single curvature (straight element) blades and includes helical 
inducers as well as the cambered M-l fuel pump (Ref. 12) , The dashed lines 
represent the tandem row designs having double curvature blading in the vane 
exit regions. These dashed lines were utilized to establish losses for the 
stage loading considered in the work split analysis. Ideal head coefficients 
for 0.29 to 0.5 were considered. The lower value is a typical one for helical 
inducers while the higher value represents a tandem row design. 

Conventional axial-flow compressor blade correlations 
of the total pressure loss parameter with diffusion factor D (Ref. 7) were 
used for the high-speed rotor analysis. The loss parameter is defined as 
follows : 

a) cos &2 
~2o 

to - loss coefficient as defined previously for 
inducer blading 

$2 _ blade exit relative flow angle (measured from 
axial direction), degrees 

a - blade solidity (chord length/blade space) 

Existing values (Ref. 15, Figure 203) for streamlines other than those near 
the rotor tip were increased by 23% to account for the tip loss effects shown 
for the outer streamlines. 

The hydraulic efficiency of the hydraulic turbine was 
estimated using a fixed total head loss coefficient (based upon exit mean 
radius relative velocity) of 0.177. In addition, a leakage factor of 0.975 
(ratio of power actually obtained to power obtained from velocity diagram 
calculations) was utilized. Experience with the design (Ref. 16) and testing 
(Ref. 17) of turbines having a moderate degree of reaction indicates that 
these are reasonable values. 

b. Interstage Power/Speed Matching Characteristics 

To evaluate the stage matching characteristics of the 
system, it was assumed that the exit relative velocity of each blade row was 
constant regardless of the operating flow coefficient. 

Writing the expression for the low-speed inducer/turbine 
energy balance at design or off-design flow coefficient assuming zero inlet 
whirl to the inducer stage as follows: 
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where : 


K U . 2 I N. 

U i Cu i + (0/0 d ) r R = TLFF | 

U - tangential rotor velocity at mean diameter, ft/sec 
Cu - fluid exit tangential velocity, ft/ sec 
N - rotor speed, rpm 

R - ratio of actual high-speed rotor/ low-speed rotor 

rotational speed to design point high-speed rotor/ 
low-speed rotor rotational speed 

K - factor to account for low-speed shaft parasitic 
losses (mechanical and shroud or film bearing 
fluid friction) , ft/sec 

0 - flow coefficient (C /U) 

m 

Cm - fluid axial velocity 

TLFF - turbine leakage flow factor 

Subscripts : 

D - refers to design point value 

1 - refers to low-speed inducer 

R - refers to high-speed rotor 

T - refers to hydraulic turbine 


D R C uR - V^i 


To simplify calculations, the second term on the left side of the equation, 
which accounts for non-pumping losses, is assumed to be independent of the 
flow coefficient but proportional to the speed squared. Fluid and mechanical 
friction effects do not follow such a speed relationship; however, the 
assumption follows the usual dynamic similarity assumptions and is sufficiently 
accurate for analyses in the region of design operating speed. 

, 2 
Dividing the above expression by to obtain a 

relationship in terms of ideal head coefficients, based upon inducer wheel 

velocity: 
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Writing expressions for the individual blade row ideal 
head coefficient-flow coefficient relationships as follows: 




Substituting the above three relationships into the 
ideal head coefficient expression: 



At the design point (0/0 ) = R = 1 and the above 

D R 

expression can be simplified and solved for N^/IL to obtain: 
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The above expression can be used together with the definition of R 


V«i - <YVd e 

to eliminate Nr/N^ from the general energy balance equation to obtain a 
quadratic expression in R, the normalized speed ratio parameter. The equation 
then can be solved for R to obtain a general expression in terms of blade row 
geometry, design point parameters, and normalized high-speed rotor flow 
coefficient : 
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The relationships derived above define general design and off-design high- 
speed shaft /low- speed shaft rotational speed ratio values; however, the torque 
speed matching characteristics of the low-speed shaft components must be 
examined to ensure that stable (hydraulically-locked) operation will be 
obtained. 


At any given operating point, it is necessary that the 
constant-speed hydraulic turbine torque coefficient-flow coefficient curve 
intersects the corresponding inducer curve at a slope which is greater than 
that of the inducer curve. This is illustrated on Figure No. 3(a). If the 
flow coefficient is displaced from the torque-matched values (resulting from 
such possible causes as flow disturbances or varying hydraulic turbine torque 
requirements resulting from rotor inertia during speed transients) in a 
decreasing direction, the speed will decrease because the inducer torque 
requirement will exceed the turbine torque and the unit will tend to return 
to the original torque and flow coefficient values. If the flow coefficient 
is displaced in an increasing direction, the speed will tend to increase 
because the turbine torque exceeds the inducer torque requirement; again, 
tending to return the unit to the original flow and torque coefficients. 

Figure No. 3(b) illustrates the opposite situation 
wherein the hydraulic turbine curve slope is less (more negative) that that 
of the inducer. In this case, if the unit is operated below the matched 
torque-flow coefficient value, the low-speed shaft will tend to increase 
because the turbine torque will exceed the inducer requirement and the 
flow coefficient will diverge from the matched value, resulting in a "run- 
away" speed condition. If the flow coefficient is greater than the matched 
value, the inducer torque requirement will exceed the turbine torque, which 
results in a decrease in speed and divergence in an increasing direction of 
the flow coefficient. Actually, steady-state operation of the inducer/turbine 
unit at the condition shown on Figure No. 3(b) would not be possible. The 
unit would either "run away" to an undesirable high speed with a corresponding 
low flow coefficient or it would not "bootstrap" and would operate at low or 
zero speed with correspondingly high flow coefficients. 

To define conditions and limits for stable operation in 
terms of blade geometry, design parameters, and operating flow coefficient, 
the head coefficient expression previously derived can be differentiated with 
respect to the normalized flow coefficient, (0/0d)r» holding the speed ratio, 
R, fixed; and then solved for (Nr/Ni)d, making use of the relationship between 
(Nr/N^d, (Nr/Nj.) and R: 
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Figure 3. - Low-Speed Shaft Torque Matching. 






This expression defines the critical value of speed ratio for stability. To 
ensure stable operation, the design speed ratio must be less than this value. 
Analysis indicated that the value of R does not vary greatly from unity for 
designs of interest over the operating flow coefficient region. In addition, 
the term K has a relatively minor influence upon the speed ratio 

R (0/0 D ) R 

value. Consequently, the stability of candidate designs was evaluated at the 
design point, although the minimum stability margin generally is obtained at 
the maximum flow coefficient* 

In estimating the required design margins to ensure 
stable operation, a calculation of the speed ratio variance was made using 
the following design parameter values: 

(Cu/U) D± = 0.425 

(Cu/U) = 0.500 

U R 

(Cu/U) DT - 0.524 

K = 0.02 

and the following estimated variance attributable to prediction uncertainties: 

Cu/U . - 10% 

K - 100% 

The calculated variance in (N R /Nj_>Q was approximately 13%. 
This value does not vary appreciably in the range of design interest and was 
used to establish stability margins. The only designs considered were those 
where the critical speed ratio for stability was at least 13% larger than the 
design value. 

The above analyses are based upon ideal head character- 
istics, assuming the blade exit relative flow angle and fluid deviation angle 
are fixed. The analysis for stability probably is conservative for conditions 


20 



where the flow coefficient is less than the design value. This is based upon 
consideration that in this region, the recirculation caused by boundary layer 
separation generally will result in higher than ideal inducer torque require- 
ments, which permit higher operating speed ratio values than those indicated 
by the ideal analysis limits. The ideal analysis is believed to be accurate 
for estimating design parameters and also for operation at above design-flow- 
coefficient conditions. 

c. Stage Cavitation 

It was necessary to examine estimated cavitation per- 
formance of the blade rows prior to selecting the stage work split to ensure 
that the contractual design requirement for cavitation-free operation of the 
high-speed rotor at 120% of the design flow coefficient would be satisfied. 

This requirement was interpreted to mean that at the 120% design flow coefficient 
point, cavitation performance of the system would be limited by the low-speed 
inducer rather than the high-speed rotor performance. It was assumed that the 
operation of the high-speed rotor at partially-cavitating (vapor present) 
conditions was permissible provided that head or torque reductions resulting 
from cavitation occur initially in the low-speed inducer as NPSH is reduced. 

The design cavitation performance requirements for the 
low-speed inducer do not directly influence the work split selection; however, 
criteria for the inducer performance at 120% of the design flow coefficient 
had to be selected to permit calculation of the available head for the high- 
speed rotor. The non-cavitating head coefficient of the helical portions of 
the inducer was estimated from available data for 6-degree to 12-degree 
inducers (Ref. 8, 10 and 11) to be 0.135 (based upon rms rotor velocities). 

This value was reduced by 10% to allow for performance loss resulting from 
cavitation. Using the loss coefficient data shown on Figure No. 2, reductions 
in NPSH caused by fluid heating were calculated for a range of flow coefficients. 
The conservative assumption that the inducer was operating at 25 ft NPSH was 
used to calculate high-speed rotor suction head. 

The total of the combined corrections for head degra- 
dation, fluid heating, and inducer suction NPSH is plotted as a function of 
inducer ideal head coefficient on Figure No. 4. This correction was applied 
to the estimated inducer head coefficient at 120% of the flow coefficient to 
obtain the available NPSH at the inducer discharge. The off-design head 
coefficient at 120% of the flow coefficient was estimated from the flat-plate 
inducer data previously referred to as well as data for the higher head 
coefficient designs (Ref. 12, 13, and 14). The estimate is shown on 
Figure No. 5 in terms of head coefficient as a function of design point 
ideal head coefficient. 

To calculate available NPSH for the high-speed rotor, 
an additional reduction of low-speed inducer discharge NPSH was made to 
account for interstage flow passage losses and fluid heating. A NPSH reduction 
equivalent to 10% of the inducer exit absolute velocity head was made. 
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DELTA ACTUAL HEAD COEFFICIENT A X 



IDEAL HEAD COEFFICIENT 


Figure 4. - Calculated Head Losses, Sum of Inlet NPSH, Inlet Cavitation, 
and Fluid Heating Effect. 
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Figure 5 


Off-Design Head Coefficient for Inducer Rotors 


The high-speed rotor inlet relative flow angle and flow 
coefficient are important design parameters because they establish the high- 
speed rotor inlet blade angle and blade height requirements. In view of the 
assumption that the high-speed rotor and hydraulic turbine would have identical 
annulus dimensions, the above parameters establish the annulus for both blade 
rows. 


To obtain the high-speed rotor inlet design parameters, 
criteria for the cavitation performance as a function of inlet flow conditions 
must be established. The available data for inducers and axial pumps concern 
the zero pre-whirl condition only. These data must be corrected and extra- 
polated for the current application where some pre-whirl exists. Available 
test data were compared upon the basis of equivalent cavitation number, K, 
(based upon inlet relative velocity) : 

Hj - H v - cj/2g 

K " Wj2 /2g 

- inlet total head , ft 

- vapor head, ft 

- inlet absolute velocity, ft/sec 

- inlet relative velocity, ft/sec 

2 

g - acceleration caused by gravity, ft/ (sec) 
and inlet relative flow angle (measured from axial direction): 

„ —1 Sml . 

6 - cos — — , degrees 


C ^i - inlet axial velocity, ft/sec 

Three methods of correlation, based upon no fluid thermodynamic head effects 
were evaluated for this purpose. In the first method, blade surface velocity 
correlations proposed by Lieblein (Ref. 7) were used to estimate blade low- 
pressure surface, minimum pressures, and corresponding cavitation number as a 
function of diffusion factor, D, and inlet relative flow angle. This 
correlation was originally developed to correlate losses with blade boundary 
layer thickness. It cannot be expected to give accurate estimates of actual 
low-pressure surface velocities; however, it is useful to obtain gross 
indications of surface velocities for empirical data correlations. The 
cavitation number can be written: 
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W,,,„ r - blade maximum surface velocity and 

MAX 

Wmax/WI can be expressed in terms of diffusion factor and discharge/inlet 
relative velocity ratio: 


W, 


MAX 

W, 


= D + 


!i 

w. 


+ 0.1 


D - diffusion factor 


W 2 - exit relative velocity, ft/sec 

The W 2 /W 1 ratio can be expressed in terms of diffusion factor, inlet relative 
flow angle, and blade solidity. The following is a comparative listing of 
experimental cavitation test data for the incipient head-breakdown point: 


Reference 

Inlet 

Relative Flow 
Angle 6 1 
(degrees) 

Diffusion 
Factor D 

© 

f w \ 

MAX ! 

TEST 

d 

( W MAxl 

RATIO 

V 

1 

[ W 1 

\ /THEORY 

18 

70.1 

0.370 

1.082 

1.234 

0.876 

19 

71.9 

0.504 

1.118 

1.261 

0.888 

20 

63.5 

0.631 

1.175 

1.326 

0.886 


The parameters are based upon the mean (rms) radius station; however, the 
diffusion factor was calculated using the ideal, integrated design head rise 
rather than the local value. An average value of test /theory velocity ratio 
of 0.883 was used to obtain data correlation curves. 


A second performance correlation which can be used is 
that of Gongwer (pages 29-40, Ref, 21): 


NPSH 

Cm 2 /2g 


1.4 + 0.085 (1 + ) 

0 
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Gongwer found that this relationship correlated the complete head breakdown 
point for centrifugal pumps with inlet relative flow angles (based upon tip 
geometry) of 7-degrees to 11-degrees. However, Wood and his associates 
(pages 929-940, Ref. 22) found that this correlation could be used to 
estimate the incipient head loss point for mixed-flow impellers with inlet 
tip relative flow angles of approximately 16-degrees. 

A third correlation, which correlates the incipient 
head loss point for helical inducers with tip inlet relative flow angles of 
3.5-degrees to 7.5-degrees, is the Brumfield (Ref. 23) criteria: 

NPSH = 

Cm^/2g 

The three methods of correlation are plotted on 
Figure No. 6. The correlation based upon the Lieblein surface velocity 
analysis has the advantage of taking into account blade fluid loading. However, 
the results of the work split analysis indicate that inlet relative flow 
angles in excess of 75-degrees are required to satisfy design specifications. 

In this region, blading of high solidity generally is used. A mean line 
solidity of 1.8 was selected for both the preliminary mean line analysis 
and the final design. This is an intermediate value between the airfoil type 
blade and the high-solidity helical, non-airfoil type blade; therefore, 
existing cavitation data for both were utilized. The solidity value of 1.8 
was used to calculate the dashed-line curves shown on Figure No. 6. Because 
of the lack of data for high solidity blading and high inlet relative flow 
angles, an envelope established by using the Brumfield and Gongwer correlations 
was used to estimate the high-speed rotor inlet geometry and flow conditions. 
These curves follow the trends established by the lower relative flow angle, 
low solidity (1.1 to 1.3) air-foil blade data correlations which are based 
upon surface velocities. Analyses of the relative inlet flow angles to the 
high-speed rotor at 120% of the design flow coefficient indicate that the 
value obtained did not vary greatly from the design point value. The curves 
shown on Figure No. 6 were used to establish the design inlet relative flow 
angle. The cavitation number, K, was calculated for the 120% design point 
because this is the point of minimum cavitation margin. 

Estimates of hydraulic turbine cavitation margins 
indicated that large margins existed for designs compatible with the low-speed 
inducer and high-speed rotor designs. Calculated design point cavitation 
numbers, based upon Lieblein's surface velocity correlations, are on the order 
of 0.3 for the designs of interest (using inlet head and exit relative 
velocities) . Operating cavitation numbers based upon estimates of head 
delivered by the high-speed rotor are on the order of 4.0. This allows more 
than an adequate margin for design and off-design operation. 
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CAVITATION NUMBER 








d. Work Split and Annulus Geometry Selection 

The low-speed inducer rotational speeds of 14,100 rpm 
and 18,800 rpm correspond to suction specific speeds of 30,000 and 40,000, 
respectively, for the design NPSH of 105 ft (25 ft NPSH) plus 80 ft of LH 2 
thermodynamic suppressionhead) , Calculations indicated that the 14,100 rpm 
designs would have unsatisfactory torque/speed stability characteristics; 
consequently, the 18,800 rpm speed was selected. 

Values of high-speed rotor ideal head rise from 1500 ft 
to a maximum value determined from the low-speed shaft inducer/ hydraulic 
turbine torque/speed stability limits were evaluated. This portion of the 
study included ideal head rise values that yield less than the contractual 
minimum over-all actual head rise values. A study of the lower head rise 
value was mutually agreed upon between NASA and Aerojet-General when pre- 
liminary calculations indicated a more desirable design might be achieved by 
decreasing the minimum to a lower value. Partial results of the study are 
summarized on Figure No. 7. The following parameters were plotted against 
the delivered over-all head rise for all three blade rows (non-cavitating) : 

Rp/Rj “ Radii ratio 

- High speed rotor and hydraulic turbine mean line 
(rms) radius 

- Inducer exit mean line (rms) radius 

- High rotor inlet fluid relative velocity 
(measured from axial direction) , degrees 

Efficiency - Over-all efficiency for three blade rows, % 

The calculations performed were for a low-speed shaft 
speed (Nj) of 18,800 rpm and a high-speed shaft /low-speed shaft ratio (SR) 
of 2.37. Figure No. 7 illustrates the results for three values of inducer 
head coefficient. Different values of both low-speed inducer exit mean line 
radius (R 21 ) and high-speed rotor relative velocity retardation factor, 
(W^/W^)^, were evaluated. 

The calculations indicate that the higher head rise 
designs have higher predicted efficiencies; however, the lower head rise 
designs have lower high-speed rotor inlet relative flow angles, which provide 
improved high-speed rotor cavitating and non-cavitating performance predict- 
ability. The lower head rise designs could have lower discharge annulus 
radii (corresponding to lower Rr/Rj values) and annulus dimensions, which 
more nearly correspond to downstream pumping elements. 

The efficiency predictions are based upon the previously 
discussed loss coefficients. The high-speed rotor inlet relative flow angles 
are based upon the limit curves shown on Figure No. 6. The Rr/Rj values are 
calculated from the high-speed rotor ideal head rise and the (W^/W^) value. 


R I 

8 
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Figure 7. - One-Dimensional Design Parameters vs Delivered Head 
R ot = 1.1R, t (W2/W1) = 0.70. 



While lower design head rise values would provide more 
conventional high-speed rotor blade geometries, it was believed that a feasible 
design could be achieved for the originally-selected design non-cavitating 
head rise of 5671 ft. The calculations indicated that maximum inducer exit 
radius and inducer ideal head coefficient as well as minimum (W£/Wi)r values 
yield the most desirable designs. The primary criterion was the magnitude 
of the high-speed rotor inlet relative flow angle. A (W2 /Wi)r value of 0.65 
was selected as the lowest yalue which provides reasonable blade velocity 
diffusion and diffusion factors. A value of inducer mean discharge radius 
10% larger than the inlet radius was selected. The inlet radius is obtained 
from cavitation performance considerations (see Section III,B,1) . Values of 
exit radius, which are significantly larger than this would result in inducer 
exit tip radii that are larger than the inlet tip radius. This was undesirable 
from the mechanical design arrangment aspect as well as for matching the 
inducer exit annulus to the high-speed rotor inlet. An inducer ideal head 
coefficient 


AH I 


where: AH^ - Inducer ideal head rise, ft 

j - Inducer exit mean (rms) wheel velocity, ft /sec 

2 

g - Acceleration caused by gravity, ft/sec 

of 0.425 was selected. This is the maximum value for which there is design 
experience with single blade rows. 

Initial calculations for the study were based upon the 
0.02 value of the low-speed inducer shroud and bearing friction coefficient. 
The final calculations were based upon final estimates for the shroud and 
bearing friction in LH 2 . The values for water will be somewhat higher. The 
coefficient is defined to be consistent with the equation format shown in 
Section III, A, 2, b. 


The selected design values are listed on Tables IV and 
V for LH^ hydraulic design and water test, respectively. 

B. DETAILED HYDRODYNAMIC DESIGN 

Hydrodynamic design is summarized on Figure No. 8, which presents 
both the geometric parameter and the velocity triangle for all three of the 
blade rows. 
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DESIGN SPECIFICATION SHEET 
HYDROGEN HYDRAULIC DESIGN 


ro 

to 

o 


05 

80 

O 

o 

o 

o 

IN 

O 

O 

o 

O 










• 

• 

• 

* 

■ 

* 

• 

* 

• 

• 

• 

« 

• 

« 

• 










O 


o 

05 


<0 

o 

IN 

o 

o 

rH 

o 

o 


o 










03 


o 

03 

H 1 

H 

80 


o 

80 

00 

to 

CO 

o 

o 












to 

to 

iH 

03 

03 


05 

05 


CO 

IN 

o 











1 


•* 




IN 



to 




to 

t*. 

























T3 
























ft 

05 
























as 

o 


i 





















o 

M 

ft 

2 

0 

o 

to 

o 

o 

















♦H 

ft 

ft 

ft * 

,# 


« 

4 

• 

to 

to 

o 

80 

IN 

00 

O 

o 

o 

o 

80 

O 

O 

o 

o 

s 

D 

M 

1-3 IN 

IN 

o 

o 

to 

05 

• 

4 

• 

* 

* 

* 

• 

* 

• 

• 

* 

* 

• 

« 

* 

ft 


S 

0- 


03 

o 

05 

03 

o 


o 

80 

to 

•H 

00 

IN 

o 

o 

03 

to 

05 

80 

o 

3 

ft 

ft 

to 



00 


to 

0 1 


o 

rH 

03 

05 

IN 


o 

o 

00 

CO 

rH 

03 

o 

>> 

Q 

ft 



1 






00 

r— ( 



o 


05 

03 


03 

to 


o 

T3 


H 




00 



1 


r» 




to 



to 






O 

ft 





pH 




to 

to 

o 

05 

IN 

03 

o 

o 

o 

o 

05 

o 

O 

O 

O 


* 

* 

4 

• 

4 

4 

* 

* 

4 

4 

4 

• 

• 

* 

• 

P 

o 


o 

IN 

to 


o 

IN 

o 

o 

rH 

o 

o 

80 

O 

P 

03 


o 

00 

03 

CO 

o 


o 

CO 

CO 

IN 

c0 

03 

o 


rf 


80 

rH 


rH 

o 


05 

00 



IN 


co 

O 

1 


•*» 




to 









00 


O 

05 

* CO ft 

ft\ g ’H 

• ft ft ft to 

^ HO 3 ft 


■ft ft 'ft OQ 01 (fi P ft ft 51 

o H 3 ft ft ft P H feO M 


TO 

X P 

* P 
ft ft I 

* • p ft 
ft ft P M 





0 



P 



ft 





ft 


s 



o 

ft 



ft 




3 



ft 



o 





05 


ft 



rH 

p 



P 




d 

05 


P 








ft 

ft 

D 



ft 

o 



0 

u 



CO 

U 


•H 



TJ 

^ — \ 





05 

H 


u 


H 

xj 


H 

3 



to 

ft 


> 



O 

ft 




3 

05 



05 

P 


05 



P 



05 

CO 

05 

ft 



CO 

ft 




O 

ft 




ft 

P 

05 

a 

P 


ft 

>s 


Jh 

CO 

CO 

O 



ft 

05 




•H 

CO 



P 

bD 

05 

ft 

3 

05 

cp 

U 

P 


ft 

o 

*H 

w 



ft 

ft 




H 




ft 

•H 

rH 

CO 

05 

rH 

ft 

05 

•H 



b 

ft 









O 

O 




05 

3 


•H 

3 

s 

ft 

OQ 


<u 

ft 


0) 




bQ 


05 

05 

3 

•H 



P 


H 

p 

o 

M 

ft 

£ 

3 


feO 


05 

OQ 



ft 

C 




C0 




P 



p 

•H 


ft 

0> 

05 


Sh 

3 

3 

n 



<15 

*H 


O 

O 


*H 


CO 

ft 

CO 

CO 

ft 

P 

CO 

H 

Q 


ft 

o 

ft 

ft 



05 



ft 

ft 

05 

O 


ft 

ft 

ft 

ft 

ft 

P 

ft 




ft 

*H 

CO 


£ 


ft 

ft 


05 

05 

> 

05 


a 

CO 

a 

a 

CO 

ft 

a 

P 

P 

05 

o 

ft 

CO 


O 


CO 



CO 

CO 

•H 

ft 









3 

3 

CD 

CO 

O 

o 

ft 

rH 



•H 

O 

u 

U 

«P 

CO 









ft 

ft 

ft 

*H 

3 

£ 

05 

ft 


O 

> 

Cl 

o 

O 

■H 










rH 

rH 

CO 

Q 

CO 

ft 

ft 


P 

•H 

ft 

05 

ft 

ft 

QQ 

3 









rH 

pH 






P 

*H 


Q 

•H 



o 

O 









05 

CD 

P 

rH 

pH 

pH 

rH 

ft 

O 

•H 


O 

ft 

-P 

ft 

♦H 









ft 

ft 

P 

ft 

ft 

ft 

ft 

&0 

ft 

O 



•H 

<H 


•P 









O 

O 

ft 

P 

P 


P 

•H 

ft 

<15 



3 

ft 

-P 

05 









U 

U 

ft 

O 

O 

o 

o 

<15 

ft 

ft 


«H 

rH 

ft 

05 

3 









ft 

ft 

CO 

H 

H 

H 

H 

=£ 

o 

CO 



ft 

CO 

S25 

C0 










.31 



DESIGN SPECIFICATION SHEET, WATER OPERATION 
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a. Geometric Parameters 



Low Speed 
Inducer 

High Speed 
Rotor 

Hydraulic 

Turbine 

Inlet Tip Dia 

7.72 

5.94 

5.95 

Inlet Hub Dia 

3.09 

4.75 

4.75 

Discharge Tip Dia 

7.40 

5.94 

5.95 

Discharge Hub Dia 

5.40 

4.75 

4.75 

Discharge Mean Dia 

6.48 

5.38 

5.38 

Number Blades 

3+3 

13 

40 

Max Blade Thk - Tip 

0.10 

0.123 

0.187 

Max Blade Thk - Hub 

0.25 

0.318 

0.187 

Discharge Flow Coefficient 

0.147 

0.156 

0.373 

Discharge Ideal Head 

Coefficient 

0.425 

0.259 

0.743 

Solidity at Mean 

2.8 

1.8 

1.43 

Leading Edge Radius at 
Mean 

0.015 

0.023 

- 

Trailing Edge Radius at 

0.013 

0.011 

0.015 


Mean 


b. One-Dimensional Velocity Triangles at Mean Radius 


U - Wheel speed 
W - Relative velocity 

C - Absolute velocity 

all units in ft/sec 


Subscripts 

1 — inlet u - tangential direction 

2 - outlet m - axial direction 


LOW SPEED c 

INDUCER 1 




U X =U 2 = 1050 


Cu x =261 Cu 2 «536 



Figure 8. - Summary of Hydrodynamic Design 
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1. Low-Speed Inducer 

The low-speed inducer speed and head’ rise were established 
as a result of the stage work split discussed in Section III, A, 2. In general, 
the maximum speed of the low-speed inducer will be limited by the suction 
performance. Higher inducer speeds mean greater cavitation margin for the 
high-speed rotor and less cavitation margin for the low-speed inducer while 
with lower inducer speeds, the reverse is true. The suction eye was sized for 
a 40,000 suction specific speed with the inlet fluid velocity head equal to 
the design NPSH of 25 ft. Then, the speed was calculated from the suction 
specific speed, flow, and the effective NPSH (the sum of the design or apparent 
NPSH plus the thermodynamic NPSH). A thermodynamic NPSH of 80 ft at 36.7°R 
was scaled from existing 37°R data (Ref. 24) and was found to be consistent 
with previous Aerojet-General experience (Ref. 6). 

The discharge mean radius was limited to a value that is 10% 
greater than the inlet mean radius. The discharge annulus was obtained from 
the blade element design consideration (viz., relative velocity ratios). 

The inlet blade angle of 6.3-degrees at the tip was obtained 
from criteria for optimum suction performance (Ref. 25) . The blade angle 
from tip to hub was established with r tan 8 being constant, which is the 
condition of a flat-plate or helix inducer. The back portion of the inducer 
was made of twisted airfoil sections. The flow was analyzed at the end of 
the flat-plate portion and again, at the discharge of the airfoil section. In 
both cases, it was assumed that the flow satisfied simple radial equilibrium. 
The flow leaving the flat-plate portion of the inducer was assumed to follow 
the blade angle without deviation. The losses were obtained from data for 
similar types of inducers (Ref. 6). 

Then, the rear portion of the inducer was analyzed, assuming 
the following conditions: 

a. The flow entering the airfoil portion was the same as 
that leaving the flat-plate portion. 

b. The mass averaged input head would be the difference 
between the over-all inducer input head (as determined by the work split) and 
the flat-plate input head, 

c. The losses would be considerably less than found in 
normal airfoil data because there would be no entrance loss. 

The over-all inducer losses would be consistent with those 
for similar inducers (see Figure No. 2) 
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The discharge axial and tangential velocity as well as the 
total and static head are shown on Figure No. 9 in relationship to passage 
height. 

With the flow conditions defined, blade stagger and camber 
angles could be established which would satisfy this flow. The inlet flow 
conditions were assumed to be the same as the flat-plate discharge flow 
condition; therefore, camber and stagger angles could be determined from the 
relative flow angles and some deviation angle distribution. Figure No. 10 
shows the inducer blade exit flow deviation angle correlation, while 
Figure No. 11 shows the smoothed setting and turning angles. 

The inducer was laid out so that the three flat-plate vanes 
blended into three of the airfoil blades. This required that the airfoil 
sections be stacked on the center-of-gravity in the radial plane but shifted 
forward (from hub to tip) ia the axial plane. The other three blades were 
considered as nearly true airfoils and would serve as partial blades. 

2. High-Speed Rotor 

The high-speed rotor supplements the work done by the inducer 
and provides the input power to the low-speed shaft. The work split analysis 
established the two major blade design criteria; the cavitation coefficient at 
the blade inlet and the blade loading resulting from the required head rise. 
Existing NASA methods (Ref. 15) were applied in the flow analysis, and was 
used to calculate the blade element performance. The input head to the high- 
speed rotor was adjusted until the mass averaged input head agreed with the 
one-dimensional input head. 

Figure No. 12 shows both the inlet and outlet high-speed 
rotor blade performance parameter. To obtain the inlet flow conditions, it 
was assumed that conservation of angular fluid momentum at any given stream- 
line existed between the low-speed inducer discharge and the high-speed rotor 
inlet. With the angular or tangential velocity known, an axial velocity 
distribution could be found which satisfied both radial equilibrium and 
continuity. This results in an axial velocity decrease from hub to tip. 

The discharge axial velocity distribution was kept nearly 
equal to the inlet axial velocity distribution so that the streamline curva-_ 
ture through the blade row could be neglected. The blade loss coefficient (to) 
was obtained from data for similar types of blades (Ref. 13). 

An existing method (Ref. 15) was used to obtain the blade 
stagger and turning angles from the fluid relative flow angles. Figures 
No. 13 and No. 14 are extensions of two figures in the referenced work 
(Figures No. 137 and No. 138, Ref. 15) and extrapolate the data into the high 
fluid inlet angle region (70-degrees to 90-degrees) . The incidence angle was 
calculated from the following: 
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Figure 9. - Low-Speed Inducer Blade Element Data (Sheet 1 of 2) 
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Figure 9. - Low-Speed Inducer Blade Element Data (Sheet 2 of 2) 








Figure 10 


Inducer Blade Exit Flow Deviation Angle Correlation 
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Figure 11. - 


Low-Speed Inducer, Stagger and Camber Angle Distribution 
for Double-Curved Portion of Blades 
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Figure 12. - High-Speed Rotor Blade Element Data (Sheet 1 of 2) 
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Figure 12. - High-Speed Rotor Blade Element Data (Sheet 2 of 2) 
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Figure 14. - Extrapolation of Slope Factor/Inlet Fluid Angle Curves 



where: 


i = i + n 4> 
o 

i is incidence at zero camber (Figure No. 13) 
o 

n slope factor (Figure No. 14) 

<|> blade camber angle 

The deviation angle was obtained from data for similar blade 
rows (Ref. 26) and is shown on Figure No. 15. With both the incidence and 
deviation angles known at each radial station, the blade stagger and camber 
angles could be calculated. Figure No. 16 shows stagger, camber, fluid, blade, 
incidence, thickness, and deviation angles at each radial station for the C-4 
blading selected. 


3. Hydraulic Turbine 

The hydrodynamic design of the hydraulic turbine which drives 
the low-speed inducer through the outer shroud is discussed in this section. 
This turbine is located directly downstream of the high-speed rotor and 
utilizes the full flow of the tandem inducer set. The design problem is to 
ascertain the blade geometry and the fluid exit conditions for a prescribed 
exit fluid angle and shaft power. The resulting design is a constant section, 
untwisted blade with a blunt leading edge. 


The following are the design criteria: 


a. The design is weighted to favor of both simplicity and 
economy in construction at the expense of hydraulic 
efficiency. 

b. Blade cavitation is avoided at design point operation. 

c. Blade cavitation at off-design conditions is not 
investigated. 

Both simplicity and economy in construction are achieved 
through the design of a constant section untwisted blade. A prime considera- 
tion in designing the profile shape is the wide variation of inlet fluid angle 
along the height of the blade. There is a 40-degree relative fluid angle 
variation between the hub and tip radii (see Figure No. 17). Consequently, 
a blunt leading edge was used to accommodate this wide variation with a 
minimum incidence loss. This type of blunt rotor profile was utilized for 
the design of the oxidizer turbopump turbine of the M-l engine (Refs- 16 and 
17). At a relatively low inlet Mach number, the design point performance of 
blunt profiles is equal or superior to sharp blading edge profiles while the 
off-design performance or operation under incidence is superior. 
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Figure 15. - Blade Exit Flow Deviation Correlation Used for High-Speed 
Rotor Blade 
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Figure 16. - High-Speed Rotor Blade Element Data (Sheet 1 of 2) 
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Figure 16. - High-Speed Rotor Blade Element Data (Sheet 2 of 2) 
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Figure 17. - Hydraulic Turbine Inlet Relative Flow Angle Distribution 
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A one -dimensional analysis was performed to ascertain the 
design point mean line velocity triangles and blade losses. The high-speed 
rotor mean discharge conditions fix the hydraulic turbine inlet velocity 
triangle. The exit velocity triangle is fixed by the shaft power requirement 
and the mass flow. A system that was developed for gas turbine design 
(Ref. 16) was utilized for the loss estimates. The blade Reynolds Number and 
Mach Number are within the limitations of the system. A constant discharge 
blade angle was desirable for simplicity as well as economy in fabrication. 
Therefore, the exit through-flow velocity was adjusted to produce a constant 
discharge fluid angle. 

4. Discharge Housing 

The discharge housing is a nonoptimum hydraulic design. In 
an actual pump application, the turbine would discharge into the inlet of the 
main high-speed pump, either centrifugal or axial. The housing in this appli- 
cation was designed to turn the flow radially in the shortest axial length to 
obtain minimum high-speed rotor shaft length. 

The guide vanes were placed to give the housing sufficient 
strength as well as to provide a passage, through which the rear hydrostatic 
bearing flow could be passed. The mean camber line angle of these vanes was 
established to match the turbine discharge flow assuring free vortex flow in 
a constant height channel. The leading edge radii were deliberately made 
large and the blade row solidity was kept low to assure that the vanes would 
not stall or have large shifts in performance at off-design flow. 

C. MECHANICAL LAYOUT 

1. Over- All Design 

Figure No. 1 shows sectional views of the major pump compo- 
nents. The drive unit, which is shown in phantom, is a slightly modified 
NERVA technology turbopump turbine/bearing housing assembly. This tester was 
driven by an electric motor for the steady-state testing and by a turbine for 
the transient testing. This turbine is the first-stage wheel of the NERVA 
assembly and gaseous nitrogen is used to power it. For liquid hydrogen testing, 
the turbine can be powered using either gaseous nitrogen or gaseous hydrogen. 

Most of the components are fabricated from 6061 Al, except for 
the high-speed rotor (7079 Al) and the rear instrumentation cluster (CRES 416). 
The 416 stainless steel was selected for the instrumentation cluster to permit 
the hardware to be used for liquid hydrogen operation with rolling element 
bearings. The 6061 Al was selected because both plate and rod c'an be obtained 
in various sizes and it has good machining as well as welding characteristics. 
The low-speed inducer is the highest stressed 6061 component and has an 
estimated fatigue safety factor of 1.59 in liquid hydrogen at the design speed. 
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Only the high-speed rotor was stressed higher than the maximum 6061 safety 
limit would allow; therefore, 7079 A1 was used. Operation at cryogenic 
temperatures will be possible because at the large radii, all nonrotating 
components are of the same material and will therefore contract together. 

The pump has a centerbody at the inlet supported by inlet 
vanes . This centerbody serves two purposes ; it allows pressure measurements 
to be taken between the high-speed rotor and the low-speed inducer and provides 
bearing support for the back-up rolling element bearings . 

A set of three total and one static pressure measurements was 
located at the high-speed rotor inlet, the high-speed rotor discharge, and the 
hydraulic turbine discharge. The total pressure was measured with Kiel probes 
at three radial stations while the static pressure was measured with a drilled 
tap at the hub shroud. The measurements were used to determine blade row 
performance during the steady-state testing. 

High-speed rotor rotational speed and torque were measured 
on the input shaft (tester unit) and the low-speed rotor rotational speed was 
measured at the shroud just ahead of the front hydrostatic bearing. 

2. Prony Brake 

The Prony Brake was designed to measure torque during steady- 
state testing. Essentially, the two rotors remain at the same speed ratio 
during steady state, even at off-design conditions. Utilizing the Prony Brake 
to apply an additional load upon the hydraulic turbine permitted the effect of 
blade row interaction at different speed ratios to be obtained. The Prony 
Brake also was used to calibrate the turbine when the inducer had been 
removed. The braking action was accomplished by means of a caliper action 
on a rotating steel disc, which was an integral part of the hydrostatic axial 
thrust ring. After the tests requiring the use of the Prony Brake were con- 
cluded, the steel ring was replaced with a 6061 A1 ring. The torque force was 
translated into a pressure reading by means of a piston-cylinder arrangement. 
The piston end attached to the brake and the cylinder end of the housing. 

3. Hydrostatic Bearings 

Hydrostatic thrust bearings and journal bearings are used to 
support the low-speed assembly (i.e. , low-speed inducer, hydraulic turbine, 
and connecting shroud). To verify the concept, a hydrostatic bearing para- 
metric analysis was conducted to select fluid film bearings that would provide 
optimum minimum flow rate and maximum load carrying performance while satis- 
fying the basic system geometry. The preliminary geometry of the inducer, the 
turbine, and the shroud assembly was designed. Estimates of the axial and 
radial forces, weights, and inertias were made so that the study would cover 
a useful range of expected operation. 


50 



The results of the parametric study showed that hydrostatic 
bearings could be used in this application. Some additional analysis was 
made for the hydrogen testing design to size the orifice diameter and the 
pocket recess depth. 

The hydrostatic bearings selected for water testing were 
analyzed with the exact design dimension. The performance of the journal 
bearings is shown on Figures No. 18 and No. 19 while the performance of the 
axial thrust bearings is shown on Figure No. 20. The calculated axial thrust 
at the design condition is 2500 lb toward the pump suction. The maximum load 
capacity toward a suction of 5000 lb is sufficient because the maximum axial 
thrust at off-design for this type of assembly normally is less than 1.5 of the 
design point thrust. Additional thrust capability can be realized by increasing 
the supply pressure. During testing, the pocket pressure could be monitored and 
the supply pressure increased whenever the ratio of pocket to supply pressure 
exceeded the design value. 

4. Rolling Element Bearings 

Provisions were made to use rolling element bearings as a 
back-up design. For the water test, these bearings were radial ball, grease- 
packed, 5200 steel bearings, with nonmetallic separators. The front bearing 
rode on the centerbody supporting the inducer end. The back bearing rode on 
the aft instrumentation cluster in the same place as the labyrinth. Both are 
shown on Figure No. 1. 


Again, it should be noted that the concept of a full-flow, 
hydraulic-turbine- driven inducer can be better designed when hydrostatic 
bearings are used on the shroud. This eliminates the inlet support vanes and 
centerbody and does not limit the inducer- turbine speed because of high 
bearing DN. 


5. Stress and Critical Speed 

The three major components (low-speed inducer assembly, high- 
speed rotor, and high-speed shaft assembly) were stress analyzed at both 
water test speed and hydrogen test speeds. The low-speed inducer assembly has 
a fatigue factor of safety of 1.59. The highest stress levels were at the 
blade root at the inducer inlet. The high-speed rotor disc has a burst speed 
of approximately 86,000 rpm, which is well above the design speed of 44,500 rpm. 
The blades show a factor of safety of approximately 7 at a total average stress 
of 6.35 ksi and an alternating stress of 0.75 ksi. The high-speed shaft 
assembly has a lateral vibration critical speed of 36,000 rpm to 42,000 rpm 
depending upon the effective shaft stiffness. If the drive turbine (NERVA 
Technology Turbine Wheel P/N 294332) were redesigned for a smaller, lighter 
wheel, the drive assembly would be capable of operating at the design speed. 
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NOTES: 


1. NUMBER OF POCKETS - 8 

2. WATER LUBRICATED 

3. SUPPLY PRESSURE - 315. PS I 

4. EXIT PRESSURE - 115. PS I 

5. TOTAL CLEARANCE - H T = .008 

6. ORIFICE DIAMETER - D 0 = .080 



Figure 20. - Double-Acting Thrust Bearing Performance 



The low-speed shroud deflection at design speed operation is 
not uniform over the width of the forward hydrostatic bearing. The change in 
deflection is negligible under the rear bearing. For liquid hydrogen opera- 
tion, the front hydrostatic bearing will either be machined conically for a 
consistent clearance on each side, or machined straight with some compromise 
in clearance and flow. This will compensate for an estimated coning type 
deflection under the bearing of 0.001-in. Deflection of the thrust bearing 
is negligible. 

D. OFF-DESIGN PERFOBMANCE 

Predicted performance curves were prepared that describe the non- 
cavitating and cavitating performance of each blade row and the discharge 
housing. These curves, along with the methods used to obtain the off-design 
data, were published in the interim report for Tasks I, II, and III (Ref. 27) 
The tabular data from these curves was used for the computer model and the 
parametric studies subsequently discussed in Sections IV and V of this report 
Similar nomenclature and format were used for each blade row and the housing 
to simplify the handling of program input and calculations. The performance 
curves obtained from the test data are included in Section VIII, C of this 
report along with the design point predictions, as applicable. 
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IV. 


TASK II— FORMULATION AND DEMONSTRATION OF COMPUTER SIMULATION 


The formulation of a digital computer model for the dynamic analysis of 
an inducer pump with a full-flow, full- admission hydraulic turbine is presented 
in this section. Two systems were modeled in conjunction with the inducer; 
a closed-loop system for transient testing with water, which was structured 
after the D-3A test facility and an open-loop system which can represent either 
a chemical rocket engine feed system or a cryogenic pump test facility 
(i.e., Test Stand H-6) . 

A. DESCRIPTION OF COMPUTER MODEL 


The general approach in modeling the inducer was to describe the 
performance of the blade rows by normalized performance curves, which were 
obtained from separate computer programs for blade and passage design. For 
example, to study a change in pitch angle, the first step is to generate 
revised performance maps. Then, the effect of the revised blade performance 
upon the dynamics would be studied using this program. 

The performance curves portraying the stages of the inducer can be 
mono-variant, bi-variant, or tri-variant functions 

\ - f( V 

Yj - g(X r Zj > 

\ - h( V V V 

The fixed tables of input values are interpolated. The work balance on the 
shafts, the calculation of boundary conditions, the energy (heat) balance, 
and the calculation of minor losses complete the inducer model. 


The system model is composed from a finite difference, method of 
characteristics, solution to waterhammer theory. Waterhammer theory is used 
to describe all hydraulic components. For right traveling waves 


JLsl + 


w 

gA 


W 


x,t 


C+ 


while for left traveling waves 

P „ V W . 

x > fc _ = c _ 

P gA p 

where C+ and C- are the waterhammer constants or characteristics. These 
equations are solved along with non-linear friction and appropriate boundary 
conditions (the inducer model being one) to describe the system (Ref. 28). 
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The pressure-weight flow rate form of the "characteristics" as a 
finite difference solution to the quasi- linear hyperbolic partial differential 
continuity and dynamical motion equations is convenient for direct comparison 
with chemical rocket engine test data. 

Appendices A and B contain the detailed assumptions, derivations, 
and equations used to portray the inducer and feed system. The techniques 
applied and assumptions made are included in the interim report (Ref. 27). 
Appendix C contains the digital program listings. The program details (i.e. , 
flow charts, nomenclature, test case, input format and output format) are 
available in Ref. 28. 

B. TASK II — CONCLUSIONS 

The following significant conclusions are drawn from the formulation 
and simulation of the computer model. 

1. The characteristic curve approach to the description of the 
dynamic performance of a multi-stage inducer pump is workable and efficient. 

Basic transients can be executed in three minutes using the IBM 360-65 system. 

2. A large quantity of data input is required by the characteristic 
curve approach. This necessitates that careful attention be paid to curve maps, 
their size, limits, and number of points input, as well as the handling of 
maximums and minimums. Emphasis must be placed upon extreme off-design char- 
acteristics because these describe the critical part of the transient where small 
errors result in large variations from the true condition. Also, careful, 
methodical verification of program input data must be made to ensure the absence 
of error. 


3. The minimum characteristic curves that are essential to evalu- 
ation of the design should be used. Other minor correctors should be set to 
zero fields. 


4. Small computing intervals are essential to the finite difference 
techniques used. 

5. The waterhammer description of the closed- loop system provides 
an analysis of the oscillatory stability characteristics of the system as well 
as the over-all "inertia" characteristics of the suction and discharge lines. 
Therefore, the program should prove valuable as an analytical tool for the 
synthesis of test facilities for dynamic pump testing as well as an analytical 
tool for existing systems. 

6. The waterhammer description of the open-loop system is 
particularly valuable for studies with hydrogen. The compressibility effects 
of rapid pressure changes in the hydrogen suction and discharge lines signifi- 
cantly affect the pump dynamics. For example, when rapid bootstrapping of the 
low-speed shaft occurs during the advance engine transient, the flow rate into 
the discharge line exceeds the outflow by 10%. 
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V. 


TASK III — PARAMETRIC STUDIES 


The inducer dynamics were investigated parametrically in two systems. 

The first system was the Test Stand D-3A test facility, \diere the demonstration 
unit was tested with water, and which is a closed-loop system. The second 
system was an open-loop system, and an advanced cryogenic engine feed system 
configuration as well as the NERVA pump test facility. Test Stand H-6, were 
investigated. Liquid hydrogen was the fluid stuided. 

A. DEMONSTRATION UNIT IN A CLOSED-LOOP SYSTEM 

The parametric study for the closed-loop system was used to determine 
the significance of different parameters upon the start and shutdown transients. 
The general closed-loop configuration is presented in Appendix D. The effect 
of each of these parameters were described in the interim report (Ref. 27), 
which was used as a guide in designing the pumping systems and the test set-up. 

The 14 different cases shown on Table VI are those showing the 
greatest change from the nominal start and shutdown condition. In each case, 
the significant parameters were compared with like parameters for the nominal 
case. 


The following conclusions are based upon the results from the 
parametric study of the closed-loop water system. 


1. The unit, as designed, will operate satisfactorily over the 
range of specified test conditions (i.e. , duration, throttling, off-design 
flow, and speed). However, because of induced oscillations in the closed- loop 
system, testing at a NPSH of less than 6.1 ft may not be feasible for most 

3 sec duration start-up tests. 

2. The induced waterhammer oscillations predicted by the model 
were experienced previously in the Test Stand D-3A test facility at similar 
frequencies. Therefore, it appears that this is a valid prediction. However, 
the actual amplitude and frequency must be verified by testing because critical 
variability could exist in the purity (i.e., lack of inert gas) of the water. 

3. The model has demonstrated that significant changes in ampli- 
tude of the induced oscillations can be obtained by utilizing larger discharge 
line sizes or increasing the inert gas content in the system. These system 
changes could be applied to control the induced oscillations during testing. 

B. DEMONSTRATION UNIT IN AN OPEN-LOOP SYSTEM 

The two,- open-loop configurations studied were the NERVA pump test 
facility (Test Stand H-6) and an advanced cryogenic engine hydrogen feed system. 
The general open- loop configuration is presented in Appendix D while the 
geometry and pressure schedules for both systems are listed on Table VII. 
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TABLE VII 

GEOMETRY AND PRESSURE SCHEDULES 
TEST STAND H-6 FACILITY 


W P i Delta PL th I.D. V <0 


Component 

lb/sec 

psia 

psi 

ft 

in. 

in. 

ft/ sec 

Tank 

29.2 

50 

0 

- 

- 

- 

- 

Suction Line 

29.2 

50 

-24 

106 

.148 

8.125 

3461 

Pump 

29.2 

26 

+57 

- 

- 

- 

- 

Discharge Line 

29.2 

83 

-6 

25 

.5 

7.625 

3487 

Valve 

29.2 

77 

-49 

- 

- 

- 

- 

Dump Line 

29.2 

28 

-4 

195 

.322 

7.625 

3481 

Outlet 

29.2 

24 

-1 

- 

- 

- 

- 



ADVANCED ENGINE 





Tank 

47 

70 

0 

- 

- 

- 

- 

Suction Line 

47 

70 

-4.00 

28 

.1 

8 

3374 

Pump* 

47 

66 

+5704 

- 

- 

- 

- 

Discharge Line 

47 

5770 

-10 

2 

.3 

3 

3464 

Valve 

47 

5760 

-90 

- 

- 

- 

- 

Injector Feed Line 

47 

5650 

-1000 

3.5 

.3 

3 

3475 

Injector 


4650 

-4627 

- 

- 

- 

- 

* Includes advanced 

cryogenic main 

stage pump 





W - Weight flow rate 

P^ - Inlet pressure to component 

Delta P - Pressure rise or drop through component 

L - Length of component 

I.D. - Inside diameter of component 

V w - Acoustic velocity through component 


60 



The detailed results from the studies of the two conf igurations , 
the effect of blade Reynolds number correction, and the examination of transient 
characteristics of the hydrodynamically similar systems were included in the 
interim report (Ref. 27). 

The following conclusions are based upon the study of the NERVA 
hydrogen pump test facility transients and advanced cryogenic engine transients. 

1. The demonstration unit will bootstrap satisfactorily for both 
a short duration (3 sec) and a long duration (9 sec) transient if tested in 
the Test Stand H-6 facility with hydrogen. 

2. No sustained induced pressure oscillations are present in 
these transients. 

3. The maximum temperature rise through the induced stages, as 
calculated using liquid hydrogen thermodynamic properties, was a 2.5°R rise 
above a suction temperature of 40°R. 

4. A rocket engine main- stage pump could be operated with the low- 
speed inducer at a minimum NPSP of 5 psi at the selected transient speed 
build-up. It is reasonable to conclude that the NPSP specified in the design 
document could be achieved by contouring the hydraulic resistance in relationship 
to time by the use of a control valve. 

5. The demonstration unit can be made to bootstrap with an 
advanced centrifugal hydrogen main-stage pump utilizing a tank-head start. 

6. For a satisfactory tank-head start with the low-speed inducer, 
precise control of the pump flow by contouring the flow control valve is required. 
The study indicates a narrow window of flow control is needed for satisfactory 
operation with a tank-head start. 
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VI. 


TASK IV— DETAILED DESIGN 


The detailed design of each test unit component was made from the master 
layout, Drawing No. 1154157. This layout included all linear dimensions, 
materials, specifications, tolerance stack-ups, run-outs, concentricity, and 
other necessary information. 

A. TOP ASSEMBLY 


The top assembly, Drawing No. 1154740, is made up of two sub- 
assemblies; the Inducer Pump Kit and the Drive Assembly. All applicable 
drawings were made in accordance with specification MIL-D-1000 (1 March 1965) . 
The top assembly consists of 12 different build-ups, which include water or 
liquid hydrogen testing, electric rotor or turbine drive, and several combina- 
tions of low-speed inducer-Prony Brake-hydraulic turbine configurations. 

B. INDUCER PUMP KIT 


The inducer pump kit, Drawing No. 1154782, consists of 141 parts, 

96 (approximately 70%) of which were standard MS, AN, NAS, or commercial items. 

As a result, the over-all drafting and fabrication costs were kept minimal. 

C. DRIVE ASSEMBLY 

The drive assembly, Drawing No. 1154352, was obtained from the NERVA 
program with some modifications. These included a new rotor shaft, grease-packed 
bearings to replace the 440C liquid hydrogen bearings, and a stub shaft that 
mounted in place of the turbine wheel. This shaft served as the means for con- 
version to the electric motor drive. The first-stage turbine wheel, driven 
by gaseous nitrogen, was used for transient testing only. 
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VII. TASK V— FABRICATION AND ASSEMBLY 


The fabrication and assembly of the pumping system were accomplished in 
the following manner. 

A. VENDOR SELECTION AND COST 

Vendor selection was based upon the application of the standard. 
Government- approved, Aerojet-General procedures wherein qualification and cost 
are the determining factors for an award. 

In this program, minimum over- all cost as well as minimal vendor 
coordination were prime objectives. Therefore, a procurement award was made 
to one supplier for the fabrication of all the pumping system components as 
well as the assembly of the first configuration to be tested. In this way, 
it was possible to readily verify the relationships between each part in the 
assembly. Also, where part modifications were found to be necessary during 
the assembly sequence, the appropriate changes were accomplished using the 
same machine and set-ups. 

B. DESIGN MODIFICATION DURING FABRICATION 

Although the test unit was designed to permit all of the components 
to be fabricated from plate, sheet, or bar stock, the vendor indicated a prefer- 
ence for using ring forgings in almost every case. This provided an over-all 
cost savings by reducing the amount of welding and pre-machining needed. A 
limited number of other changes also were made (i.e., materials for bolts and 
inserts as well as bolt head configuration) because difficulty was experienced 
in obtaining the originally specified item. 

Both the low-speed inducer and high-speed rotor blades were made by 
utilizing tracing patterns. The inducer pattern was made full size and the 
high-speed rotor pattern was made ten times full size. The hydraulic turbine 
blades, which were a constant section, were directly machined. All of the 
fabrication was completed in approximately three months. 

C. DESCRIPTION OF MAJOR HARDWARE 

The low-speed inducer is shown on Figure No. 21 encased in the 
rotating shroud. The combined weight of this inducer and shroud is 8.06 lb 
and it has an inertia of 0.254 lb-in. -sec^. The concentric rub marks which 
are evident on the photograph were caused when the high-speed shaft bearing 
failed (see Section VIII, B for detailed discussion). 

Figure No. 22 is a view of the high-speed rotor showing the blade 
leading edges as well as machining markings. This rotor weighs 1.94 lb and has 
an inertia of 0.0178 lb-in. -sec2. 
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The braze joint between the shroud and blade tips of the hydraulic 
turbine is shown on Figure No. 23. Again, concentric rub marks also are 
evident. The combined weight of the hydraulic turbine and shroud is 3.06 lb and 
inertia is 0.088 lb-in. -sec2. 

Both the forward hydrostatic thrust bearing and the forward 
hydrostatic journal bearing are shown on Figure No. 24. The removable orifices 
for both of these bearings as well as the aft bearings, which are not shown, 
were all of the same diameter. Removable orifices were used so that if necessary 
additional load capacity could be obtained by installing orifices of different 
diameter. 


Figure No. 25 shows the instrumentation centerbody and the probe 
plate, which mounts between the low-speed inducer and the high-speed rotor. 

This is a non-rotating assembly. The total Kiel probes and the wall static 
probe at the hub can be seen in this view. The plug, which can be identified by 
two staking marks, can be replaced with a diaphragm type high-frequency pressure 
transducer. The taper that mates with the inlet guide vane housing is located 
at the end opposite the probe plate. 

Figure No. 26 is a view of the Prony Brake assembly with the thrust 
ring mounted. The brake assembly remains stationary in the pump housing while 
the thrust ring is both attached to and rotates with the low-speed inducer/ 
hydraulic turbine shroud. The actuation pressure is supplied to six pistons 
through 3/ 16- in. tubing. 
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Figure 24. - Forward Hydrostatic Thrust and Journal Bearing. 



Figure 25. - Instrumentation Centerbody 
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VIII. TASK VI— SYSTEM TESTING AND MATHEMATICAL MODEL REFINEMENT 


This section describes the water-loop testing and the refinement of the 
mathematical model. The testing phase with the pumping system served to pro- 
vide performance data for each blade row under steady-state, cavitating, and 
non- cavita ting conditions as well as to simulate a main-stage pump start and 
shutdown transient. The performance data then was used to refine the mathe- 
matical model wherein predicted performance data had been used previously (see 
Section III,D). 

A. WATER-FLOW-LOOP FACILITY 

Test Bay D-3A (see Figure No. 27) at Aerojet-General’s Sacramento 
Facility was used for the water-flow- loop testing. Figure No. 28 is a flow 
schematic of this test facility. Some facility modification was accomplished 
to permit concurrent operation of Test Bays D-3A and D-3B. This was necessary 
because an auxiliary flow supply was needed for the hydrostatic bearings and 
to actuate the Prony Brake. A second-stage Titan oxidizer pump was installed 
into Bay D-3B and it delivered an auxiliary flow of approximately 70 gpm at 
600 psi. This flow was filtered and the pressure reduced so that the flow to 
the bearings in the test unit was maintained at 250 psi above the suction 
pressure. The bearing discharge flow was returned to suction upstream of the 
facility flow straightening section. 

Two methods were used to drive the test unit. The first method 
was used for all steady-state testing. It consisted of a 1000 hp synchronous- 
speed motor delivering power through a variable-speed eddy current clutch and 
two, speed-increasing gearboxes. The maximum output speed of this motor drive 
was 8500 rpm, which was the speed that most testing was conducted. The second 
method consisted of replacing the motor drive with a GN£ driven turbine. This 
turbine was a first-stage NERVA wheel that was operated at a maximum speed of 
9000 rpm. The turbine speed was regulated by a closed- loop servo-controlled, 
hydraulically-actuated GN 2 flow valve. 

B. TEST PLAN AND SCHEDULE 

The test schedule was designed to assure maximum test facility 
utilization by this program as well as two other concurrent contractual 
efforts; the Hubless Inducer Program (Contract NAS 3-11218) and the Thrust 
Balance Program (Contract NAS 3-7978) . This resulted in the testing of one 
program being conducted while hardware was being rebuilt to a different con- 
figuration in another program. 

1. Type and Number of Tests 

The pumping system was tested 11 times in three basic types 
of tests; head-flow non- cavitating, cavitating, and transient. There were 
five head-flow non- cavitating tests, four cavitating tests, and two transient 
tests. Table VIII shows the type of test planned for each build-up along with 
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TABLE VIII 


TYPE TEST FOR EACH BUILD-UP 


Build-Up Type 

Test No. P/N Test Input N 

1,2,3 1154740-2 H-Q 8500 

6800 

8500 

6800 

8500 

6800 

6800 

4, 5 1154740-2 CAV 8500 

8500 

8500 

6800 

6 1154740-3 H-Q 6800 

6800 

6800 

6800 

7 1154740-1 H-Q 3400 

5100 

6800 

8500 


Fixed 


Low- Speed 

Valve 

Variable Valve 

Maximum 

PS = 85 

Q = 935, 

747, 840, 


psia 

935, 1030, 1121, 

Maximum 


maximum 
Q = 595, 

747, 897 

2520 


Q = 745, 

840, 935, 

2010 


1030, 1120, maximum 
Q = 590, 745, maxi- 

1440 


mum 

Q = 745, 

840, 935, 

1050 


1030, 1120, maxi- 
mum 

Q = 590, 745, 840, 

0 


935 

Q = 590, 

745, 840, 

Maximum 

Q =935 

935, 1030, 1120, 
maximum 

PS = 50, 30, 20, 15 

2520 

1440 

0 

Q =935 
Q =935 
Q =745 

10, 8, e 

-2, -4, 
psig or 
loss 

> 

> , 4 , 2 , 0 , 
-6, -8, -10 
until 10% 

f 

Maximum 

PS =85 

Q = 450, 

, 595, 747, 

2010 

1050 

0 

Maximum 

PS =85 

897, 1050, maximum 

i 

Q = 150, 225, 300, 



375, 410, 450, 
maximum 

Q = 225, 335, 450, 
560, 615, 675, 
maximum 

Q = 450, 600, 750, 
825 , 900 , maximum 
Q = 655, 745, 935, 
1030, 1120, maximum 
1120, 1030, 935, 
745, 655 
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TABLE VIII (CONTD) 

TYPE TEST FOR EACH BUILD-UP 



Build-Up 

Type 



Fixed 



Test No. 

P/N 

Test 

Input N 

Low-Speed 

Valve 

Variable Valve 

8-9 

1154740-1 

CAV 

6800 

Maximum 

Q=600 

PS = 50, 

30, 20, 15, 







10, 8, 6 

, 4 , 2 , 0 , - 







-4, -6, 

-8, -10 psig 







or until 

10% loss 




6800 


Q=750 






6800 


Q=900 






8500 


Q=750 






8500 


Q=935 






8500 


Q=1120 

> 

f 

10-11 

1154740-4 

Trans 

0 to } 
8500 

No 

PS to 

Q = 0 to 

750 - 3 sec 




Control 

be de- 

trans 



termined 


8500 

to 


Q = 750 to 0 - 

3 

see 

0 



trans 



0 to 



Q * 0 to 935 - 

3 

see 

8500 



trans 



8500 

to 


Q = 935 to 0 - 

3 

sec 

0 



trans 



0 to 
8500 

(1) 

No 

Control 

PS=Min Q = 0 to 935 - 
trans 

6 

sec 

8500 

to 


Q = 935 to 0 - 

6 

sec 

0 



trans 



0 to 



Q - 0 to 935 - 

3 

sec 

8500 



trans 



8500 

to 


Q = 935 to 550 

- 

2 

5500 



sec trans 



5500 

to 


Q = 550 to 935 

- 

1.5 

8500 



sec trans 



0 to 



Q = 0 to 935 - 

6 

sec 

8500 



trans 



8500 

to 


Q - 935 to 550 

- 

4 

5500 



sec trans 



5500 

to 


Q = 550 to 950 

- 

3 

8500 



sec trans 
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the data to be recorded. The actual amount of data recorded was somewhat 
greater and under slightly different conditions. Most of the testing was 
conducted at the maximum motor drive speed of 8500 rpm; however, some data 
points were taken at lower speeds to observe anomalies in speed scaling and 
to obtain data at zero flow conditions. 

All non-cavitating test suction pressure was set at 85 psia 
to assure that no cavitation existed in any of the blade rows. The flow was 
varied by a butterfly type valve located in the pump discharge line. A maxi- 
mum 125% of the design flow could be obtained. This upper limit was reached 
when the over-all pump discharge pressure rise equaled the flow loop pressure 
drop with the butterfly valve wide-open. 

In the cavitating tests, the flow-to-speed ratio was main- 
tained at some predetermined value and the suction pressure lowered until the 
pump had either reached 10% to 20% head loss or the delivered discharge pres- 
sure was insufficient to flow the loop with the butterfly valve wide-open. 

The start transient tests were conducted by manually adjusting 
the speed control system. In this system, a variable resistance potentiometer 
set the speed at which the drive turbine would run. The shut-down transient 
tests were conducted by closing the GN 2 supply valve and overriding the speed 
control system. 


The hydrostatic bearing supply flow always was set at the 
operating condition before the test pump was rotated. This was done to keep 
the bearings from rubbing although a workable system could be designed to 
supply the bearings with a pressure consistent with a main-stage pump discharge 
pressure. However, a system of this type would require the bearings to rub 
during the initial start-up. 

2. Instrumentation and Recorded Data 


Table IX is a listing of the test instrumentation and includes 
the parameter symbol, description, range, and location along with the method(s) 
applied to record the data. 

Both the inlet and discharge pressures were measured using a 
manifold with four taps located 90-degrees apart. All other pressure measure- 
ments were made by means of a single tap located internally in the pump. Three 
total Kiel probes and one static wall tap were located at the discharge of each 
blade row. The Kiel probes were situated at 20%, 50%, and 80% of the passage 
height. They were set at an angle that was based upon the design prediction. 
These probes indicate a nearly-true pressure over a flow angle change of 
+ 20-degrees. The static wall taps were 0.030-in. in diameter and located at 
the hub. The set of probes at the low-speed inducer discharge actually was 
located closer to the high-speed rotor inlet (see Figure No. 1). Consideration 
was given to locating probes at the inducer discharge, but was discarded because 
of the probe length and design complexity. Drilled passages connected both the 
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TABLE IX 


INSTRUMENTATION LIST 


PARAMETER 

DESCRIPTION 
A- BAY 

FULL-FLOW INDUCER 

RANGE 

PS-IF,2F, 

3F,4F 

Static Suction Pressure 

-12-100 

psig 

PSI HF-2 
or 

PSHF MS-1 

Static Suction Pressure 

-12-100 

psig 

TS 

Suction Temperature 

0-200°F 

PD 

Static Discharge Press. 

0-200 

psig 

PD 

Static Discharge Press. 

0-200 

psig 

TD 

Discharge Temperature 

0-200°F 

PBS 

Hydrostatic Bearing 
Supply Pressure 

0-300 

psig 

TBS 

Hydrostatic Bearing 
Supply Temperature 

0-200°F 

QBSF and 
QBSR 

Hydrostatic Bearing 

Supply Flow (Front and Rear) 

0-100gpm 

PTBP-1/ 

Pressure of Front Thrust 

0-100 

PTBP-2 

PTBP-3 

Bearing Pocket (3 Loca- 
tions ) 

psig 

Q 

Flow Rate 

0-1500 

gpm 

NHS 

High Shaft Speed 

0-10,000 

rpm 

NLS 

Low Shaft Speed 

0-4000 

rpm 

PPBS 

Pressure Prony Brake 
Supply 

0-300 

psi 

IDS-1 

Inducer Disch. at Hub 
Static Pressure 

0-100 

P§ig 

PIDT-1 

Inducer Disch. Total 
Pr. at R = 2 *506 -in. 

0-100 

psig 

PIDT-2 

Inducer Disch. Total 
Pr. at R = 2.691-in. 

0-100 

psig 

PIDT-3 

Inducer Disch. Total 
Pr. at R = 2.867-in. 

0-100 

psig 



RECORDER 


< 

as 

3 

H 

a 

% 

3 

j as 


TD I*} 

H 

J*. w 

w ft 


10 o 

K 

o 

0, 

^3 

o a 

LOCATION 

W <tj 

> O 

h as 

CO O 

M 

Q 

as g-i 

Upstream St. Vanes 
Manifold 4 Taps 


X 

X 



Upstream St. Vanes 
Close Couples Protocon 
or M.S. Transducer 




X 

X 

Upstream St. Vanes 



X 



One Pipe Dia. Down- 
stream from Discharge 
Flange Manifold 3 Taps 


X 

X 



Discharge Line, Photo- 
con or M.S. Transducer 




X 

X 

One Pipe Dia. Down- 
stream from Discharge 
Flange 

X 


X 



Downstream of Flow 
Control Valve 

X 


X 



Downstream of Flow 
Control Valve 

X 





Downstream of Hydro- 
static Bearing^ Flow 
Control Valve 


X 

X 



Adapter Housing Front 


X 

X 



Adapter Housing Front 



X 



Adapter Housing Front 



X 



In. Main Flow Loop Line 


X 

X 

X 

X 

Jack Shaft or Test Assy 
Hsg. (Turbine Drive) 


X 

X 

X 

X 

Adapter Housing 


X 

X 

X 

X 

Upstream of Prony Bk. 
Actuation Valve 

X 


X 



Inlet Housing 



X 



Inlet Housing 



X 



Inlet Housing 


X 

X 


X 

Inlet Housing 



X 






TABLE IX (cont.) 






RECORDER 

PARAMETER 

DESCRIPTION 

RANGE 

LOCATION 

VISUAL 

GAGE 

STRIP 

CHART 

DIGITAL 

H. FREQ. 
TAPE 

OSCILLO- 

GRAPH 

PIDHF-1 

Inducer Discharge Pr. 
High Frequency 

0-100 

psig 

Inlet Housing Photocon 
or M.S. 




X 

X 

PIDHF-2 

Inducer Discharge Pr. 
High Frequency 

0-100 

psig 

Inlet Housing Photocon 
or M.S. 




X 


PRDS-1 

Rotor Discharge Static 
Pr. at Hub 

0-300 

psig 

Adapter Housing Rear 



X 



PRDT-1 

Rotor Discharge, Total 
Pr. at R - 2.506-in. 

0-300 

psig 

Adapter Housing Rear 



x 



PRDT-2 

Rotor Discharge Total 
Pr. at R - 2. 691 -in. 

0-300 

psig 

Adapter Housing Rear 


X 

x 



PRDT-3 

Rotor Discharge Total 
Pr. at R = 2. 864-in. 

0-300 

psig 

Adapter Housing Rear 



X 



PTDS-1 

Turbine Discharge Static 
Pr. at Hub 

0-200 

psig 

Adapter Housing Rear 



X 



PTDT-1 

Turbine Discharge Total 
Pr. at R - 2. 506 -in. 

0-200 

psig 

Adapter Housing Rear 



X 



PTDT-2 

Turbine Discharge Total 
Pr. at R = 2.691-in. 

0-200 

psig 

Adapter Housing Rear 


X 

X 


X 

PTDT-3 

Turbine Discharge Total 
Pr. at R = 2.864-in. 

0-200 

psig 

Adapter Housing Rear 



X 



PACC 

Static Accumulator 
Pressure (if needed) 

-12-100 

psig 

Suction Accumulator Photo- 
con or M.S. Trans. 




X 

X 

PFCV 

Control Valve Discharge 
Pressure 

0-100 

psig 

Downstream of Main H^O 
Flow Control Valve Photo- 
con or M.S. Trans. 




X 

X 

LFCV 

Control Valve Position 

0-100% 

Main HgO Flow Control 


X 


X 

X 

PTF-2 

Torque Force Converted to 
Pr. 

0-1000 

psig 

Pump Housing 


X 

X 



T 

Rotor Torque 

0-1000 

in-lb 

Jack Shaft 

X 

X 

X 




DRIVE TURBINE (WHEN USED) 








PDTI 

Pressure Drive Turbine 
Inlet 

0-100 

psig 

Inlet GN^ Line 

X 


X 

X 

x 

TDTI 

Temp. Dr. Turbine Inlet 

-200 + 
100°F 

Inlet GN,> Line 



x 


x 

LVTI 

Turbine Valve Position 

0-100% 

Inlet GN2 Line 


X 

X 

X 

x 


B-BAY 

(AUXILIARY PUMP) 








PS 

Suction Pressure 

-12-100 

psig 


X 





PD 

Discharge Pressure 

0-500 


X 





N 

Drive Speed 

0-10,000 

rpm 



X 




Q 

Flow 

0-1000 

gpm 



X 




TBPE 

Bearing Temperature 

0-300 °F 

Pump End Roller Brg. 

X 





TBT 

Bearing Temperature 

0-300°F 

Thrust Bearing 

X 








total Kiel probe and the static wall tape sensing ends to the transducers 
located on the outside of the pump housing. The inducer discharge pressure 
readings were incorrect during the first five tests because of leakage in the 
drilled passages at the junction of the inlet guide vane and the suction spool 
housing. This area of leakage was bridged by pressing drilled dowels into the 
passages at the junction interface. 

Visual gage and strip charts were used to record critical or 
"redline" parameters as well as aids in the setup of a test. The strip charts 
also provided "quick- look 1 * information and served as a backup for the digital 
recording. 


Aerojet-General’s Automatic Data Control System (ADC) was used 
to record the digital data. Each parameter or channel is "swept" 30-times-per- 
second and the data is recorded on magnetic tapes. Three-second time sweeps 
were used during steady-state testing while continuous sweeps were utilized 
during the transient operation. The steady-state data for each function were 
averaged and statistically analyzed over the three-second time interval. This 
statistical analysis of each function included the standard deviation and the 
maximum difference from the mean. Then, the data were examined and considered 
to be acceptable if the 2-sigma deviation was less than 1% of the mean. 

The above described average values were computed using Job 
No. 3001. This computer program, which has both a printed and punch card out- 
put, changes the digital counts to engineering units. The punch card output, 
which was identical with the average values recorded on the printed output, 
then was used as input to a data reduction program. The test data were 
utilized in this data reduction program to produce normalized pump performance. 

High-frequency data was recorded for the purpose of determin- 
ing the pressure oscillation during cavitation. However, limited cavitation 
data was obtained (see Section VIII, C) and consequently, the data from only 
one cavitation test was reduced. These data did not show any significant pres- 
sure fluctuations; therefore, the remaining data were not analyzed although 
they were recorded. 

Oscillograph data were used for "quick- look" purposes to 
ascertain if the start and shut-down transients had fulfilled contractual 
requirements. Those transient graphs included in this report were obtained 
from the ADC recorded data and plotted using a digital data plotter (Cal-Comp 
Plotter) . 


3. Pump Build-Ups 

As indicated, the initial pump build-up was completed during 
the fabrication phase. All of the remaining build-ups were accomplished in or 
adjacent to the test bay. Those build-ups requiring less than a normal working 
shift to complete them were done without removing the pump from the test bay. 
When major hardware changes were necessary, the pump was removed from the test 
stand to permit the bay to be used for another testing program. 
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The turbine-end bearing of the NERVA power transmission failed 
during the first test because it was of a type other than that called out on 
the drawing. The bearings were specified under the FF-B-171 military system 
and this one had been misordered in converting to commercial parts. 

The originally specified bearings were grease-packed with non- 
metallic separators. It was found that they could not be obtained within the 
constraints of schedule requirements; therefore, bearings that were compatible 
with an air-oil mist lubrication system were substituted in the NERVA power 
transmission. These bearings were available commercially and the lubricating 
system was on-hand in the test area. This bearing system proved to be com- 
pletely successful during the remaining tests. 


When the bearing failed in the first test, some damage was 
incurred because the high-speed rotor blade tips rubbed against the low-speed 
shroud. This rub damage was removed by making a light machine cut (0.005-in.) 
on the outside diameter of the high-speed rotor. The damaged area inside the 
low-speed shroud was removed to a depth of 0.015-in. and then built back up 
with epoxy, which then was oven-cured. The shroud was then re-machined to 
provide the necessary design tip clearance between the rotor and the shroud. 
Damage to the outside diameter of the low-speed shroud was limited to some 
light rubbing contact with the inside diameters of the front and rear hydro- 
static bearings. Only the shroud required some manual attention with a crocus 
cloth. Both hydrostatic bearings functioned properly throughout the entire 
test program. 

C. TEST RESULTS 

All of the test results discussed in this section were obtained 
from the water-loop testing, which was completed between 8 October 1968 and 
23 January 1969. The three basic types of tests conducted were non-cavitating 
head-flow, cavitation, and transient. 

The test results presented for each blade row are based upon the 
mean Kiel probe reading unless otherwise noted. The blade performance is 
based primarily upon the mean Kiel probe reading because it provides a better 
indication of the blade characteristic at off-design flows. The pressure dis- 
tribution at the low-speed inducer discharge (high-speed rotor inlet) and the 
high-speed rotor discharge show considerable variation from hub to tip, with 
the tip pressure being the higher (see Figures No. 29 and No. 30). As expected, 
the discharge pressure of the turbine was nearly uniform from hub to tip even 
at off-design conditions. 

1. Steady-State Non-Cavitating and Cavitating Tests 

The steady state non-cavitating head flow tests (Tests No. 1, 

2, 3, 6, and 7) were conducted with the suction pressure sufficiently high to 
preclude any cavitation. The Prony Brake was used in Tests No. 1, 2 and 3 to 
vary the speed ratio between the high-speed and low-speed shafts. These tests 
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Figure 30. - High-Speed Rotor Passage Head Distribution 
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were for the purpose of obtaining the high-speed rotor and the hydraulic tur- 
bine performances at off-design speed ratios. 

Test No. 6 was conducted to determine the hydraulic turbine 
performance. In this test, the low-speed inducer was removed and a false shroud 
was substituted to support the low-speed rotating assembly. This testing pro- 
vided the ’key' which enabled the low-speed inducer and hydraulic turbine power 
to be separated. Without the low-speed inducer, all of the hydraulic turbine 
output power could be equated to the measure Prony Brake power plus the shroud 
friction power. This shroud friction power, which was calculated to be 60% of 
the hydraulic turbine discharge power at the design condition, is made up of 
the forward and aft hydrostatic journal bearings, the forward and aft hydro- 
static thrust bearings, and all the remaining surface on the outside diameter 
of the shroud. It was assumed that the inside diameter shroud friction would 
be less and could be included in the hydraulic turbine efficiency. This assump- 
tion appears valid because part of the shroud covers the hydraulic turbine and 
part of it covers the high-speed rotor, which would produce reverse drag. 

The hardware used in Test No. 7 more nearly duplicated that of 
an actual pumping application. The Prony Brake had been eliminated, but the 
instrumentation clusters remained between the rotating blade rows. With this 
configuration, the speed ratio was determined entirely by the hydraulic power 
match between the blade rows. The speed ratio was approximately 16% below that 
of the design value (0.426), partially as the result of the larger shroud 
friction. The design speed ratio of 0.426 had been calculated using hydrogen 
shroud losses. 


Tests No. 4, 5, 8, and 9 served to provide steady-state cavi- 
tation performance for each of the blade rows. The design requirement speci- 
fied that the over-all performance would not be limited by rotor cavitation up 
to 120% of flow. This was interpreted to mean that the low-speed inducer could 
cavitate provided that the head produced from the inducer precluded cavitation 
in the high-speed rotor. However, inducer cavitation characteristics could not 
be obtained because the high-speed rotor was the first to cavitate under all 
flow conditions. Low-speed inducer and high-speed rotor test data show three 
reasons for the high-speed rotor cavitating before the low-speed inducer. These 
are: 


- Low-speed inducer head rise lower than design 

- Speed ratio lower than design 

- High-speed rotor cavitation performance lower than design 
These reasons are discussed in the sections that follow. 

a. Low-Speed Inducer Performance 

The normalized head-flow and torque-flow of the low-speed 
inducer are shown on Figures No. 31 and No. 32, respectively. The shape of the 
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Figure 31, - Normalized Head-Flow of the Low-Speed Inducer 
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normalized head flow curve was similar to the predicted value; however, the 
low-speed inducer head at the design flow was approximately 50% of the pre- 
dicted value. This performance shift could be attributed to any or all of the 
following three reasons: 

(1) The data shown was taken from a measurement station 
located near the high-speed rotor inlet. This means that the actual discharge 
head of the low-speed inducer could have been higher and the loss between the 
low-speed inducer and the high-speed rotor were larger than predicted. As can 
be seen on Figure No. 1, there is a reversal in the passage between these two 
blade rows and most certainly, three-dimensional flow exists. 

(2) The pressure measured by the mean Kiel probe may not 
represent the mass averaged head. Figure No. 29 shows the measured head at 
three different radial stations but the flow distribution is unknown; there- 
fore, only a simple average can be taken from the three probe readings. This 
average results in a 20% increase in measured head but assumes that the flow 
distribution was constant from hub to tip. 

(3) The low-speed inducer was designed for a theoretical 
head coefficient of 0.425. The data shown on Figure No. 2 indicate that this 
is the maximum head coefficient for a single blade row. If the flow through 
the low-speed inducer was not two-dimensional as assumed in the design analysis, 
a flow shift could occur which would overload the airfoil portion of the blade. 

The normalized torque-flow curve does not include the 
low-speed shroud friction. The data shown must be added to the shroud torque 
to obtain a power balance with the hydraulic turbine. 


The suction performance of the low-speed inducer, shown 
on Figures No. 33 and No, 34, was taken from the design prediction because 
cavitation data could not be obtained for the inducer. This performance is 
needed as input to the transient model when operating in a cavitating condi- 
tion. The method used to generate these predicted characteristics were 
described in detail in the interim report (Ref. 27) . 

b. High-Speed Rotor Performance 

The normalized head-flow performance of the high-speed 
rotor is shown on Figures No, 35 and No. 36, Performance at one speed ratio, 
that test value closest to the design point, is shown on Figure No. 35 along 
with the predicted design point head. The actual delivered head was 12% higher 
than the predicted value. This higher head could be accounted for partially 
by the speed ratio being lower than the design value. The stall point was 
somewhat difficult to define because of the interaction between the blade rows. 
The shape of the performance curve was steeper than predicted, which indicates 
that the blades were "choking off" at high flows. 
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Figure 33. - Low-Speed Inducer H/NPSH vs S 
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Figure 34. - Low-Speed Inducer T/ p NPSH vs S 
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The normalized head-flow at various speed ratios is 
shown on Figure No. 36. Data points are not shown because of the cross- 
plotting necessary to obtain discrete speed ratios. These curves are some- 
what different from the design in that the effect of speed ratio was less 
significant than predicted. 

High-speed rotor efficiency in relationship to normal- 
ized flow is shown on Figure No. 37. This efficiency is approximately 5 
percentage points higher than predicted. 

Cavitation performance of the high-speed rotor is shown 
on Figure No. 38. Three different flows are shown at essentially the same 
speed ratio and two different speed ratios at the same flow. The high-speed 
rotor cavitation performance at complete breakdown and design flow was 9000 
as compared to the predicted value of 25,000. The problem here is a basic 
lack of design criteria for this type of blade. Cavitation parameters have 
not been developed for airfoil type blades with inlet fluid angles greater 
than 80-degrees. The data used in the design were extrapolated from blades 
where the testing was at inlet fluid angle of 70-degrees. Better cavitation 
performance could have been obtained if the design head-rise would have been 
less which results in lower blade loading. Improvement also would have been 
realized if the blades had been double-circular arc rather than the thicker 
C-4 type. 


c. Hydraulic Turbine Performance 


The normalized head-flow and torque-flow for the 
hydraulic turbine are shown on Figures No. 39 and No. 40, respectively, with 
speed ratio parameters. At a zero speed ratio, the values of head and torque 
divided by speed squared would be infinite; therefore, a speed ratio of 0.01 
was used instead of zero. This assumed that the measured head drop and 
measured torque would be the same at 70 rpm as at zero rpm (the speed at which 
the measurements were taken) . 

The shape of the head-flow curves are similar to those 
predicted. The difference occurring at the high flows where the predicted 
curves returned to zero head (zero head at normalized flow of 0.35). If the 
data are extrapolated to the design speed ratio of 0.425, the normalized head 
agrees within 10% of the predicted value. 

The normalized torque-flow curves had a characteristic 
that was different from the prediction although the design point agreed 
reasonably well. As was discussed in the interim report (Ref, 28) there were 
little data available to account for off-design hydraulic turbine performance. 


91 



iii»!5!isss5sss«!!»5"s"«s”""*!”ss:sssssssss:s:ss:ss5ssss5ssl 


■■ ■■■■■■■■■■■! 
■•■■■■■■■■■■■■I 
■■■■■■■■■■■■■■I 


■■>■>■»■>>>> >h 

|«MM>a»MBaal 

■■■■■■■■■■■■■■I 


RM«MRI«H| 

IHIII 


■■■■■a 

■«■■■■ 

■■■■■■ 

■■■■■■ 

■■■■■■ 

■■■■■a 

■■■■■a 

■a 

■a 

■a 

55 

aaa 

■■■ 

iiiiia 

■■■■■« 

aaaaaa 

■■■■■a 

■■■■■■ 

55 

■a 

55 


■■■■«W 

■■■■aw 


■»»a»aaaa^H 

lijiSiiSjaa 


■■■■■■ 

[■■■■■■i 

Will. 


■■■■■■ 

■■■■■■ 

■■■■■■■ 


555 : 55:25551 

ibiisiiisiipI 

■■■■*■■■■■■■ 

■■■■■■■■■■■I 

■■■■■■■■■■■■I 

■■>>>■■»>»>■ 

ssssssl 


I::::: 


55” >55 I 

■■■■■■I 


■■■■■■I 


■■■■■■■■■>■■■■■■■■■»■ 
■■■■■■■»■■■»■■■■■■■«■ 
imiiiiiiiiiMiiiiiiii 

Wmmmmwmm aaaaaaaaaaaaa 
■■■■■■■■■■■■■■■■■■■■■■a 
■■■■■■■■■■■■■■■■■■■■•■ 
>>»»»■>■»■■■■■■■■■ 
■>■■■■■■■■■■■■■■■■■■■■ 
»»»»»»»»■»■■■ 


■■■»■■■»■>■»■>■■«■■ 

»>>»>»■»»>>»■>» 

■■■■■■■■■■■■■■■■■■■■■■■i 


:::::::::::::::::::::::::::: SiHBsSSsaaSiSSiSi 


a aaiaffl iaiaaigiBjjg 


>555555555 555 > ! 



i>»»»»»»i| 
l»»>»»>»»l 
!■>>■>»»■■■■>■ 


■>>»»»■■ 

■■■■■■■■■■I 


»»»>»»■ 


555555555555555S5555S555i55555555555>5555555555»H B aa55aa 

5555555555i5555555>>«>i>>»»»»aBaBBBB»*Ba*"*"""""ii"»" 


5r^555555S55555S5S555S55SS!5SKK5K5555KS5SSS5S5SS555KSl 




55555555555555iiiiiiiiiiiiii>jjj>>>>>» 4 «>««i 


«■■■■■■■■■■■■■ 

55555555555555 

15555555555555 

5555555555555 


!■■■■■■ ■■■■■■■■■■■■■■■■'!■■■■■■*£ 

M— — ■ — — lyiyiy-— Bi 
illlllMHI 


■ ■■■ 
5555 


■■■■■■■■ 
. ■§■■■■■>» 
■ •■ ■■■■■■■■ 


»>■■■■■ 

■■■■■■»>■■» 
■■■■ ■■■■■■■■ 


5555 


I ■■■■■■ 

K ■■■■■■ 


■■■■■■■■■■■■■■■■■■■■■■■■■■55! 


■■■■■■■■■■■■■a K-'SJ SSSKS! 

Si!i! SSSB 3HliilB 

iaassassRi 


888:8»8:i88!s:::| 


■«■*■■«■■■■■> «■■■■■■ 
!**■•■- '***■*! 


i!i;KiKiB K!K;;;!;;!ai::g 83i3iS888;888l 


f «■■■■■■■■■ 


555555! 




iiiiiiniiiirnfBpiHHa 

» ■:< 5nn5r5Tjp>>>>>> 


Ba5aa5aaaaaaaB»Baariaaaaaaa«aa»*aiBagggggggEE 
■■■■■■■■■■■■■■«■■■■■■■■■■■■■■■■£■■■■■■■■■■■■ 


aaaaa55555r,55Sii5555K55KK K555555K 

ISS^SSssasssBBsa 


'SS5"5KSS5S5S^^5555S55555SS5r5S55«SK5Sw5 




■■■■■■ 


■«■■■■ 

5S555 


aii>W> M«M«>««««««M>MM>> 

■5555555555555555555555555555 

I \ !SS5SSS5SS555 55555 55555 
■»>»»»» 

i>»«»>MM 


■■■■■■•■■■■■■■■■■■■■■■■■■■■■■■■■■■■■I 
Hi>>>>» >»»>«» >«>>>>>>»»>>«■»■ 
i>>MM>>>M>M*>M»>M>5>*55!!"5»5"5E55*2”5*5*5552l 

BSIKSSKSSiSmiaSSii SSiKiKSKSKiKiSijiiKSSi l 
i $ ^^^^555555555555^^^ « * wm mm 


■■■■■■■■■■■■■■■■■< 


!■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■•■I 
■ *■■■■■■■■■■■■■■«■■■■■■■■■■■■■■■■■■■■ 
!■■■■■>■■■■■■■■■■■■■■■■■■ ■■»■■■*■■■■■■ 

Mllllllllll'llllllllllll|l||lll|lil|| 
■■■■■■■■■■■■■■■■■■■■■■ 


liiimi 

I aawia* 
| ■■■■■■! 
I >■>■■■ 
I «■■■■■ 


I nm 


■■■■■■■■■"■■■■■■■aai 

(■■■■■■■■■■■■■■■■■■■l 


555 


■■■■■■■■■■■■■■■■■■■■a 


555555555555555555555 

£55555555555555555555! 

»»»»»»»»■■■» 


•»»■ «■■■■■■■■■■■■■■■■■■■«■ 

MUM 5555555555555555555555 
* ^5555555555555555555555 
55555S5mmi m 


mmHMiiuilMIH 

■■■■■■■■■■■■■■■■■■■■■■ 

■■■■■■■■■■■■■■■■■■■■■■■I 

>»>»»»>>»»»»» 

■■■■■■■■■■■■■■■■■■■■■■ 


■■■■■■■■■■■■■■■■■■■■■■ 


5555555555555 

■■■■■■■■■■■■■ 


5555555555555 555555 
5555555555555MMI 


■■■■■■■■■■■■a 

■>■■■■■>■■■■■ 




HjHHHHjHjjlHjHlUjHHjHjOiUs::::::::::::::::;;:::;::: 


!55555555555555555aaB55555aaw«aaaa»»a« giBaagaaagggggggggggggBgggggE i 

laaaHHaaaaaaaaaaaaaBM 


SSSaEiiSS 

i»»h»»»» ■■*■■“■■*■"*“““■■■52225252225 


8888: 888 

MMiaaMil 


■■■■■■BBBBBriHta 

55555S5555555a5555555555555S5555aB>>>>>>a>aB»-“«*******»«2S“222 

■SBSS55S 55555555555555555555555555555555555555aawiaaaaaaaaaaaa«a«a 
■1225>>B>>>>>nB>>aaa ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■ 
55555555555555S55555S555555555»»>>>>>bb>>>>-->«---»»-- 


■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■“■■■■■■■■■55222552222222I 
55555S555555555555555555555aSaaaaaaaaaaaaa*aaaaa>aa»>a»a>| 
■5555555555555a555Ia55aaa5aaaaaaaaig H “***“ BB "**“ M,l liliiiilll 

>>>»»■ »■■■!»« ««■»■»■ 


>>»»»»>>■ 


( ■»»>■»>>>»»>»>■ 
»>»a>aaaaa»a»a| 


Pi*ggggggggi5“5Bgg5“l 


Q/N f gpm/rpm 


Figure 37. ^ High-Speed Rotor Efficiency vs Normalized Flow 
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Figure 38. - Cavitation Performance of the High-Speed Rotor 



Figure 39. - Normalized Head-Flow for the Hydraulic Turbine 
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Figure 40. - Normalized Torque-Flow for the Hydraulic Turbine 
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d. Over-All Performance 

Figure No. 41 shows the over-all pumping system perfor- 
mance. The system characteristics were normalized using the delivered flow 
and the head difference from suction to discharge, as well as input speed and 
torque. Again, the design point for each parameter is marked for reference. 

These data were taken from the high-speed (8500 rpm) testing and represent the 
best over-all characteristic. 

2. Transient Tests 

The transient tests (No. 10 and No, 11) were conducted to obtain 
the transient characteristics of the pumping system. All of the necessary tran- 
sients are listed on Table VIII. Test No. 10 was conducted at noncavitating 
conditions while the same transients were repeated in Test No. 11 with a lower 
inlet NPSH which resulted in a partial high-speed rotor cavitating condition. 

Table X is the nomenclature for the transient graphs, includ- 
ing the unit and description for each parameter. Eight different transients 
are provided as Figures No. 42 through No, 49, These were taken from Test 
No, 10 and are typical of the more than 44 separate transient data swings taken 
during the two tests. The following is a listing of the condition and elapsed 
time for each data point. 


Data 

Point 

High Shaft 
Speed (rpm) 

% Q/N 

Elapsed 
Time (sec) 

1 

0 to 8500 

0 to 80 

4.5 

2 

8500 to 0 

80 to 0 

6.0 

4 

8500 to 0 

100 to 0 

6.0 

5 

0 to 8500 

0 to 120 

3.8 

6 

8500 to 0 

120 to 0 

6.0 

7 

0 to 8500 

0 to 100 

4.4 

10 

8500 to 0 

100 to 0 

4.7 

13 

0 to 8500 

0 to 100 

3.1 


To minimize the number of figures, only the speed-flow characteristics of most 
of the transients are shown. The pressure data, which are not shown were con- 
sistent with the speed and flow represented on Figure No. 41 

To obtain the transient data at any particular condition, 
the input speed was held at the desired value and the flow control valve 
trimmed until the required flow was obtained. Then, the valve was left in this 
position and the input speed returned to zero. Following this, the input speed 
was manually ramped to the desired speed in the required elapsed time. This 
means that the over-all pumping system was essentially operating along a con- 
stant specific speed line. To obtain a specific speed at the end of a swing 
that was different from the start of the swing, both the input speed and the 
flow control valve were manually ramped. In most cases, these manual ramps 
resulted in the typical "S"-shaped transient. 
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Figure 41. - Over-All Pumping System Performance 


97 



TABLE X 


NOMENCLATURE FOR TRANSIENT GRAPHS 


PARAMETER 

UNITS 

DESCRIPTION 

Q/N-I 

gpm/ rpm 

Pump flow/low-speed inducer speed 

Q/N-R 

gpm/rpm 

Pump flow/high-speed rotor speed 

NI/NR 

rpm/rpm 

Low-speed inducer speed/high-speed 
rotor speed 

NI 

rpm 

Low-speed inducer speed 

NR 

rpm 

High-speed rotor speed 

Q 

gpm 

Pump flow 

PS 

psia 

Sue ti on pressure 

PD 

psia 

Discharge pressure 

PD- 1 

psia 

Low-speed inducer discharge pressure 

PD-R 

psia 

High-speed rotor discharge pressure 

PD-T 

psia 

Hydraulic turbine discharge pressure 

PP 

psia 

Hydrostatic thrust bearing pocket 
pressure (indication of axial thrust) 
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Figure 42. - Data Point 1 (Sheet 2 of 2) 
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Figure 43. - Data Point 2 (Sheet 1 of 2) 
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Figure 43. - Data Point 2 (Sheet 2 of 2) 



103 


Figure 44. - Data Point 4 (Sheet 1 of 2) 
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Figure 44. - Data Point 4 (Sheet 2 of 2) 
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Figure 45. - Data Point 
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Figure 46. - Data Point 
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Figure 47. - Data Point 
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Figure 48. - Data Point 10 
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Figure 49. -Data Point 13 


Figure No. 49, which is a fast start (3.1 sec) at the 
nominal specific speed, illustrates the pumping system characteristic. The 
high-speed rotor rotates for approximately 1 sec before the low-speed shaft 
responds. In the next second of elapsed time, the low-speed shaft has reached 
80% of the steady-state, hydraulically-locked speed ratio. During this time, 
the input speed has reached only 53% of the steady-state speed. Another way 
of stating this condition is that the two rotating shafts are 80% locked at 
slightly more than 50% of input speed. 

The condition of over-all constant specific speed occurs 
at approximately the time that the normalized flow of the rotor reaches a 
stable value. Figure No. 49 shows this to be less than 50% of the total 
transient. Up until this time, the fluid inertia in the water loop causes 
lower values of normalized flow. The flow rate increase was almost linear 
with time for most of the start transient. 

D. MODEL REFINEMENT 

Refinement of the computer transient model was the final technical 
effort in the program. The purpose of this refinement was to obtain better 
agreement between predicted transients and actual test transients. The input 
speed ranges of the transient test were not controlled (viz., the value of 
speed versus time was not pre-programmed). It was necessary to use the test 
value of speed in relationship to time as input to the computer program so that 
both the actual test transient and the computer program transient had the 
identical input speed versus time characteristic. 

1. Change in Input 

In the simulation and parametric studies previously completed 
with the computer program, predicted values of design and off-design perfor- 
mance were used for the low-speed inducer, high-speed rotor, and hydraulic 
turbine. The data for all three blade rows from the steady-state testing (see 
Section VIII, C) now were used as input to refine the computer model. These 
performance characteristics of each blade row plus the measured value of 
inertia for the low-speed shaft were the only refinements made. 

It was found that the normalized torque-flow and normalized 
head-flow of the hydraulic turbine were the two most sensitive inputs to the 
computer program. Figures No. 50 and No. 51 show these parameters on an 
expanded scale so that the complete map can be seen. As discussed previously, 
the data shown at a speed ratio of 0.01 was actually taken at zero speed ratio. 
The results remain valid because the measured torque would not be significantly 
different if the hydraulic turbine had been rotating slowly. The transient 
test data path from Data Point 1 is overlayed on Figures No. 50 and No. 51 to 
show where the hydraulic turbine operates during transient. It is interesting 
to note that the speed ratio changes from 0.01 to 0.15 in an elapsed time of 
0.5 sec while in the next 0.5 sec, it continues to increase to 0.210. Locating 
these points on the transient path line shows that the most sensitive region 
is at speed ratios of less than 0.2. 
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The computer program required a linear intrapolation between 
curves; therefore, it was important that input data be at closer intervals 
than was obtained from the test data. Consequently, it was necessary to 
extrapolate the test data to complete the other speed ratio curves. 

2. Simulations Compared to Test Results 

Figure No. 42, which contains Data Point 1, was selected as 
the comparison between the test transient and the computer model transient. 
The computer transient is indicated by a geometric symbol while the test 
transient is marked with an arable number. 

The most critical and important parameters are considered to 
be delivered flow and delivered pressure (viz., hydraulic turbine discharge 
pressure) . These two parameters are used in the design and the analysis of 
the downstream or main-stage pump. It can be seen that the computer predic- 
tion of Data Point 1 agrees very well with the test results for these 
parameters. 


The computer predicts the low speed shaft (NI) rotation to 
occur approximately 0.25 sec ahead of the actual test rotation. Most likely, 
this is caused by breakaway torque in the test unit. The computer does not 
have a logic built in for such a characteristic. Refinement to include this 
phenomenon was not considered to be significant and a value could not be deter- 
mined. Actually, inclusion of such a characteristic without a measured or 
calculated value would have been relatively meaningless because it would have 
meant that input data had to be manipulated until the computer and test data 
transients were in agreement. 
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IX. PROGRAM CONCLUSIONS AND RECOMMENDATIONS 

A. The full-flow inducer concept provides a feasible hydrogen boost 
pumping system for rocket engines, based upon the design and off-design 
analyses (including transient analyses) . These analyses also are applicable 
to other propellants with boost pumping systems of equivalent specific speed. 

B. Low-speed inducer head rise capability and low-speed inducer/ 
hydraulic turbine torque/speed stability characteristics place an upper 
limit upon the high-speed shaft/ low-speed shaft speed ratio values that can 
be used. 


C. Systems with lower design head rise and lower rotational speed than 
those adopted in this program would have more conventional blading and could 

be more easily matched to the suction geometry of main-stage pumping systems. 
This system is convenient for experimental investigations because the 
relatively large diameter exit annulus provides for back-up rolling element 
bearings and instrumentation systems; however, this annulus size is greater 
than found in axial or centrifugal main-stage pumping system inlets. 

D. The contractual requirement for cavitation-free operation of the 
high-speed rotor at 120% of the design flow coefficient is a rather severe 

one because of the relatively steep low-speed inducer head flow characteristic. 
This results in comparatively large high-speed rotor inlet relative flow angle 
(measured from the axial direction) design values which are outside the range 
of existing axial airfoil blade, cascade, compressor, and pump test data. 

Lower angles can be used if the pump cavitation performance is allowed to be 
limited by the high-speed rotor rather than the low-speed inducer cavitation 
performance at flow coefficients of less than 120% of design. 

E. Cavitating and non-cavitating test data for airfoils having vari- 
ous blade profiles and blade thicknesses as well as non-airfoil blading in the 
range of inlet relative flow angles from 75-degrees to 85-degrees (measured 
from the axial direction are needed to eliminate some of the uncertainty 
associated with the design of the inlet portions of the high-speed rotor. 
Generally, this blade row must operate with some pre-whirl and the available 
pump test data are for zero pre-whirl; therefore, a requirement also exists 
for investigating geometries designed for inlet pre-whirl. 

F. It is estimated that the inducer design for this program has the 
largest, reasonable ideal head coefficient for an axial discharge rotor of a 
single blade row design. Tandem blade row designs also can be used; however, 
they were not selected because it was expedient to retain simplicity in the 
design and fabrication of the initial demonstration unit. The tandem blade 
row design can furnish larger over-all head coefficients and the use of such 
head coefficients simplifies the high-speed rotor inlet design. Also, the 
flatter tandem row head-flow characteristics provide improved above-design- 
flow coefficient cavitation performance for the high-speed rotor. For these 
reasons, the tandem blade row should be considered for full-flow, hydraulic- 
turbine-driven inducers . 
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G. To simplify design, the high-speed rotor and hydraulic turbine use 
axial-flow, constant annulus dimension blading. Generally, systems with mixed- 
flow, variable annulus dimension blading will be required for rocket engine 
application. A mixed-flow hydraulic turbine, wherein the exit annulus is 
smaller than the inlet annulus, should be evaluated because this type of design 
more readily matches the downstream elements. 

H. For high head-rise applications, which includes this program, the 
use of more highly- loaded, high-speed rotor staging offers design advantages 
because the annulus dimensions can be made almost equal to those of the main 
stage pump suction annulus. Tandem blade rows can offer advantages for such 
applications; however, this type of design was not selected in the program 
because of the need to simplify design and fabrication. The tandem blade row 
has a disadvantage in that its added length and weight can provide high-speed 
shaft assembly critical speed problems. Compact (short axial length), highly- 
loaded staging would have application for the high-speed rotors of full-flow, 
hydraulic-turbine-driven inducer systems. 

I. Analytical investigations of the transient start-up and shutdown 
as well as the throttling characteristics of the full-flow, hydraulic-turbine- 
driven inducer system indicate that satisfactory operation was obtained. The 
low-speed shaft speed does lag behind the high-speed shaft speed during a 
start-up transient. However, for a very rapid start (<1 sec) bootstrapping of 
the low-speed shaft could be prevented by this lag. Analyses were based upon 
the type of starting experienced by a rocket engine pumping system using a 
centrifugal pump in a so-called "tank-head" start, wherein the pumping system 
is required to operate at relatively high (above design) flow coefficients 
during the initial portions of the transient. It is recommended that other 
types of operation, including those wherein the pumping system initially 
operates at low or zero flow coefficient, be considered for future investiga- 
tions. In addition, the use of other main-stage pumping systems, including 
axial flow pumps, should be considered. 

J. The test operations as well as the test results demonstrated that 
a system of this type has potential for boost pump application. Control was 
maintained easily and at no time were there any erratic or unstable operating 
conditions. However, the cavitation level of the system must be verified 
before it could be used as a boost pump. 

K. The computer transient model demonstrated acceptable performance 
prediction of an actual test start transient. Only the rotating blade row 
characteristics were changed from the original model to obtain these predictions. 
From the transients analyzed, it appeared that little if any improvement could 
be made in the prediction accuracy by additional program refinements. The form 
in which the blade row characteristics are input could be altered so as to make 
them more acceptable for the computer model, 

L. The conclusions and recommendations associated with the details of 
formulating the computer model as well as the results of the parametric analysis 
were itemized at the end of each discussion (Sections IV and V). 
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APPENDIX A 

PUMPING SYSTEM SUBROUTINES 
(BASIC EQUATIONS) 
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Appendix A 


I. NOMENCLATURE 


Parameter 

A 

B 

C 

Cp 

d 

E 

f 

g 

H 

J 

K 

Kr 

Kw 

L 

N 

P 

Q 

R 

T 

t 

t* 

th 

Vw 

w 


Description 
Flow Area 
Bulk Modulus 

Waterhammer Characteristic Constant 

Heat Capacity 

Inside Diameter 

Modulus of Elasticity 

Function of 

Gravitational Constant 
Head 

Polar Moment of Inertia 
Constant bias 

Flow Recirculation Constant 
Valve Flow parameter 
Constant Length 
Speed 

Pressure - Unless Otherwise Spec. 
Flow Rate 

Hydraulic Resistance 

Temperature 

Time 

Past Times 

Thickness of Conduit 
Sonic Velocity in Fluid 
Weight Flow 


Units 



Psi 

Ft 

Btu/lb °R 
In. 

Psi 


32.2 ft/sec^ 
Ft 

? 

Ft-lb-sec 


Lb^^-in. /sec 

Ft 

Rpm 

Psia 

Gpm 

9^9 

Sec /ft -in. 

°R 

Sec 

Sec 

In. 

Ft /sec 
Lb /sec 
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Appendix A 


Parameter 


e 


AH 

V 


Ap 

A(H/N 2 ) 
A(p/pN 2 ) 
AN/ At 
AQ/At 
AQ 

NPSH 

P 


Subscripts 

1 

2 

I 

M 

R 

T 

C 

NC 

NOM 

R 

e 

V 

A 


Description 

Error 

Total Differential Head 
Torque 

Correction to Torque 

Correction to Dimensional Head Coeff. 

Correction to Dimensional Torque Coeff. 

Shaft Acceleration 

Flow Acceleration 

Flow Rate Correction 

Net Positive Suction Head 

Weight Density 

Description 

High Speed Shaft 

Low Speed Shaft 

Low Speed Inducer 

Collector or Main Stage Pump 

High Speed Rotor 

Turbine 

Cavitation Condition 
Non-cavitation Condition 
Nominal 

Reynolds Number 
Sat Vapor Condition 
Assumed Value 


Units 


Ft 

Ft-lb 

Ft- lb 

Ft/ (rpm ) 2 

Ft-lb/lb/ft 3 (rpm ) 2 

Rpm/sec 

(Gpm) /sec 

Gpm 

Ft 

Lb /ft 3 
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Subscripts 

e 

i 


Description 

Exit 

Inlet 


r 


Recirculation 


tr 


Transient 
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Appendix A 


INDUCER SCHEMATIC AND MODEL 


P ST p SM 

P S ■> % Psi Pdx ?SR ^DR %T j t 



Inducer Rotor Turbine Collector 


BASIC ASSUMPTIONS 

1. Constant weight flow rate - no storage. 

2. Flow in positive direction only. 

3 . Rotation in positive speed direction only. 

4 . All characteristic curves dealing with head rise 
or loss are expressed as total head or enthalpy. 

5 * Nj is known as a function of time. 
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Appendix A 


II. INDUCER SUBROUTINE (J=l) 

A. ASSUMPTIONS 

1. All head heat and torque characteristics are independent 
of high-speed rotor speed, N^, 

2. No recirculating flow. 

3. A correction in stage flow and head rise exists because 
of Reynolds Number. 

4. Constant density through the stage, inlet density is 
used, Pj. 

5. Flow through the inducer is in a positive direction only 

<«j. 0) . 

B. INDUCER: ORDER OF SOLUTION 


Equation 


Boundary Conditions 

No. 

Equation 

or Limits 

1. 

Predict a Value of i , - lb/sec 

i >_ 0 

2. 

^ 448 . 831 to , = gpm 

Q = — ~Wh or 



P 


3. 


4. 


5. 


6. 
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Appendix A 


Equation 

No„ 

7, 

8 . 

9 . 

10 . 

11 . 

12 , 

13, 


Equation 

Using first past time value of N „ 

N r N * 1 

2 2 


S = 




(NPSH) 


3/4 


AH O 1 

NPSH = f 2 ( - S * N 2 *^ 

AH = K 0 (NPSH) f. 
c 2. I 

(See Curve 67-116) 


AH 


nc 


N * 
2 


f (-2 — ) 

3 v N 2 * 


AH = K, N * f, 
nc. 3 2 3 


(See Curve 67-114) 


AH V .1 

_ f (Si — Si _ ) 

AH 5 'N * ’ Q m 

nc 2 Nom 


AH 


(See Curve 67-120) 
1 


nc 


AH Re 

AH nc (1 " K 5 ah * 


nc 


AH 


AH -AH - ft 
nc c 


dl 


P S + 


AH p 
144 


psia 


Boundary Conditions 
or Limits 


0 < S < MAX 


£ 1 MAX 


0 < Q /Q. t < MAX 
— ^Nom — 


AH. 


- 1 f K, 


Re 


- 5 AH 


< + 1 


nc 
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Equation 

No. Equation 

BEFORE CONVERGENCE: BYPASS TO ROTOR 

CONVERGENCE: CONTINUE 

AP C Q 1 

14 * p (NPSH) f 6 (S ’ N 2 * ) 

(See Curve 67-117) 

Ap = K, p NPSH f , 
c 6 6 


15. 

1 

b . _ f 

Q 1 

( n 2 a) 

2 7 

P H 2 * 


(See Curve 

67-115) 


M 1 - K_ p N. 

/ k 

* 2 f 
> 7 

16. 

1 A 

p = y -Ap c 


17. 

Test Option i 

1 2 Flag 


If H 2 Flag - 0 Incompressible 

Go to Equation 19 
If H 2 Flag = + Compressible 
Go to Equation 18 
COMPRESSIBLE 

18, Compressibility effects are taken into 

account by changing density from stage 
to stage using a-* NBS thermodynamic 
subprogram.. 

The suction density to the stage is 
taken as the stage density. 


Boundary Conditions 
or Limits 


SUBROUTINE. AFTER 
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Equation 

No. 


Equation 

For Inducer Stage 

T is calculated from Heat Transfer 
Subroutine HTSR(H) 

Call PT 


f(P 0 


. T 


t-1 

JL 

En g = enthalpy - Btu/lb 
Ens * 

For all the other Stages 
Call PT 

P J = f(Pd (J-l) t _ 1> Te (J-l) t ^ 1 


^ ^dC T-] 


8 - d ^t-l, e(W) H 


For all the Stages 

The exit enthalpy if calculated 
from the torque 


AE 


773 «j Oj 


Boundary Conditions 
or Limits 


En = En + AEn 
e s u 

The exit temperature and density are 
then calculated from thermodynamic 
properties 

T = f(En , P.) 
e e d 

p f (En , P ) 

e e d 
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Equation 

No, 


19, 


Equations 

Boundary Conditions 
or Limits 

INCOMPRESSIBLE: 

rp rp I 1 

Calc. T~, P 

e 

,^ L N * 

M T30T J AH 

T •> T 
c — o 

1 — i + 7 “ 

e o C 

P 

77^ Q jPj * 778 

% * (T e> 



RETURN 
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III. HIGH-SPEED ROTOR SUBROUTINE (J=2) 


A. ASSUMPTIONS 

1, All head and torque characteristics are a function of speed 
ratio, ^/N^ 

2, Recirculation exists in the pump 

3, A correction in stage flow and head rise exists because of 
Reynolds Number. 

4, Constant density through the stage, inlet density to the 
high-speed rotor is used. 

5, Flow through the rotor is in a positive direction only, 

6, A correction to head and torque characteristics may exist 
because of high-speed shaft acceleration. (Optional) 

7, First past time value of ^ is used. 

B. HIGH-SPEED ROTOR: ORDER OF SOLUTION 


Equation 

No. 


1 . 


2 . 


3, 


Equation 

P * P. = P T (Incompressible) 

R eI (t-l) 

Q = ... ^ Q = Q = gpm 

R R(t-l) 1 

P = P 
SR dl 

Using First Past Time Value of = PH* 


Boundary Conditions 
or Limits 


4. 

q r = 

K 

r 

(Ah*) 





At 

.. Q (t-1) " Q (t-2) 

At 



5. 

* Q Re 

Q 


f 5 

1 v At’ Q„ ) 
Nom 

AQ/At <_ 

MAX 



(See 

Curve 67-146) 

q/q Nom 

MAX 
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Equation 

No, 


6, 


Equation 


Q 1 - Q [l + K x — + Q 


7. 


T s t -1 T el t _ 


t-1 


NPSH * 


P = f, (T S) . ,) 
v te b t~l 

144 (P -P ) 
sr v 


9 * 


10 , 


S - 


AH 

c 

NPSH 


vi_ 

(NPSH) 3/4 
= f, (S, 


_ai !l 

n i 5 N i 


(See Curve 67-123) 
AH = K„ NPSH t n 


11 . 


AH _1 N * 

B£ _ £ (-2_ j— 

2 3 N ’ N 

N, 1 1 


(See Curve 67-121) 


AH » K, N f 
nc 3 13 


Test Option Flag B 
If Flag B » 0 

Go tc Equation 13 
If Flag B = + 


Boundary Conditions 
or Limits 


If Q 1 < 0 Error 
Message = 140 


NPSH > ,01 


0 <: S < MAX 


0 < 4- MAX 
~ N 1~ 

N * 

0 < ~±~ < MAX 
— N, — 
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Equation 

No. 


12 . 


13. 


14. 


AH 


N 


Equation 

tr f . AN 1 M . 
-7 ■ f 4 ( TT’ V 

(Optional Curve) 


AH 


tr 


K 4 N 1 f 4 


AH 


re 


AH 


= f* C 


1 N * 1 

c _ si_) 


5 V N 1 ’ N, 9 Q„. 
nc 11 ^Nom 

(See Curve 67-127) 

AH 

AH 1 = AH (1 - K. -T7~) 
nc nc 5 AH 

nc 


Boundary Conditions 
or Limits 


AH 

•i< k 5 ^< 1 

nc 


15. 


16, 


AH 


17, 


18. 


AH - AH - AH 
nc c tr 


AH p. 


P „ + 


dR ‘sR 144 

BEFORE u)a CONVERGENCE; BYPASS TO TURBINE SUBROUTINE# AFTER 
CONVERGENCE: CONTINUE 

Au 


— f fq 

p(NPSH) r 6 N 1 9 


(See Curve 67-124) 


Ap = K, P (NPSH) f, 
c 6 6 


1 n 1 N 0 

± — s f (5_ _J) 

o 2 7 'N 9 N ’ 

P N, 1 1 


(See Curve 67-122) 


N- 


K ? p N x f ; 
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Equation 

No. 


19. 


20 , 

21 . 


Equation 
Test Option Flag B 

If 0: Go to Equation 20 

If +: Continue 


Au AN. A 

__t£ = f f __l N 

M 2 8 '-At » 1* N/ 

pN ' 1 


(See Curve VIII) 
4u tr ' K 8 p K l 2 £ 8 


M 51 ^ ' AM C ■ A “tr 


Test Option Flag 

If Flag = 0, Incompressible 

Go to Equation 23 

If Flag = 1, Compressible 

Continue 


22. COMPRESSIBLE: Calc, Te, Ene p 

e 

Methods is the same as Outlined under 
Inducer 


Boundary Conditions 
or Limits 


RETURN 


INCOMPRESSIBLE: 


T = T 
e 's 


+ 


C 

P 


Calc T , p 
e e 

r re , T 

u (30) N 1 
778 di 


AH 

778 


(T ) 


T •> T 
e — o 


RETURN 
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IV. TURBINE SUBROUTINE 

A. ASSUMPTIONS 

1. All head and torque characteristics are a function of speed 
ratio, ^/N^. 

2. A bypass flow exists in the turbine. 

3. A correction in stage flow and head rise exists because of 
Reynolds number. 

4. Constant density through the stage, inlet density to the 
turbine is used. 

5. Flow through the turbine is in a positive direction only. 

6. First past time value of N 2 , N^ = N^* is used. 

B. TURBINE: ORDER OF SOLUTION 


EQUATION BOUNDARY CONDITION 

NO. EQUATION OR LIMITS 


1 P = P T - P er ( t « X ) (Incompr . ) 


2 

3 

4 


5 


6 


Q . Q =* Q r * CPM 

H q T pT(t-l) 


p S P 

ST dR 


Using past time value of aH t = AH* 


Q r - -K t (I AH* I) 1/2 


| AH* I absolute value of AH* 
AQ/At = (Q(t-l) - Q( t-2))/ At 


AQRe 


fl ( 




At’ Q. 


) 


•NOM 


(See Curve 67-146) 

Q 1 - Q (1 + K x + Qr 


AQ/At «. MAX 

V%m - ^ 
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EQUATION 

NO. 


EQUATION 


T = T 
sit -13 erft:-l) 


Pv = f Te (1 *tt-l)' 


NPSH = 144 


(Pst - Pv) 


Using first past time value a 


10 


N 2 * 


(NPSH) 


3/4 


N * 

AHc _ f JL_ 2 v 

" 1 2 V^J XT XT ) 


(NPSH) 


V’ n i 


(See curve 67-116 for type) 


AH = K (NPSH) f 
cl l 


11 


AH „1 

- f, (2- 


N 


2 * 


»2 ‘ 3 V \ 


AH = K, N * f, 
nc 3 2 3 


(k 3 = - X.XX) 


AH n l N * .1 

re _ f /Q_ _2_ Q_ \ 

AH " 5 V N * ’ N. * Q. ' 

nc 2 1 NOM 


(See Curve 67-130) 


AH 


nc 


AH 

AH (1 - K, AU 
nc 5 AH 


re, 


nc 


i /. 


AH = AH - AH 
nc c 


BOUNDARY CONDITION 
OR LIMITS 


NPSH - .01 


0 _<_ S MAX 
0 < §- <■ MAX 

-n 2 - 

N 2 

0 < < MAX 

- Nl - 


This will be a 
negative value, 
at steady -state. 


AH 


-1 < K 


re 


- 5 AH ~ 
nc 


-b 1 
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EQUATION 

NO. 


15 


16 


17 


EQUATION 

AH 


BOUNDARY CONDITIONS 
OR LIMITS 


dt 


P + 

ST 144 


BEFORE Wa CONVERGENCE: BYPASS TO MAIN PUMP SUBROUTINE 


AFTER CONVERGENCE : CONTINUE 


Aye 


(NPSH) 


1 N * 

f c _2_\ 

f 6 N 2 *’ Nj, ' 


(See Curve 67-117 for type) 
Ape - K 6 p (NPSH) f 


PN 2 * 


_ V 

7 v N 2 *2, N x 


C 


) 


18 

19 


(See Curve 67-128) 

u 1 = K ? pN 2 * 2 f ? 

1 A 

U - M - Auc 


Test Option Flag 

If H flag « 0, incompressible go to equation 
20 1 

If flag » 1, compressible - continue 

COMPRESSIBLE: CALC. T, E , Pe 

c ne 

The equations are the same as described in 
the Inducer Subroutine except, the definition 


AEn 

u 


m T 30 N 2 
778 Q f Pj 


EXIT 
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EXIT 


A 


w 1 


BOUNDARY CONDITIONS 
OR LIMITS 


T •> T 
e — o 


HP = horse power 


s 2 -° 
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V. MAIN PUMP SUBROUTINE 


The following two options were considered: 

Option I: Main pump as collector — no rotating parts. 

Option II: Main pump as a pump. 

A. ASSUMPTIONS 

1* All head and torque characteristics are a function of speed 
ratio, N^/N^. 

2. Recirculation exists if used as a pump. 

3. A correction in stage flow and head rise exists because of 
Reynolds number. 

4. Constant density through the stage, inlet density is used. 

5 e Flow through the rotor is in a positive direction only. 

6. A correction to head and torque may exist because of high 
speed shaft acceleration. (Optional) 

B. MAIN PUMP: ORDER OF SOLUTION 


Equation 

No. 


1 . 


2 . 


Equation 


0 P M P eT 


Ct-1) 


Q = Q 


(t-1) 


M p 


Q t = GPM 


M 


(t-1) 


Boundary Conditions 
or Limits 


3 . 


P 


SM 


P 


dT 


4. 


5„ 


= T 


’t-1 


eT 


t-1 


NPSH = 144 


f Te (T eT t _ 1 ) 
(P SM “ 


NPSH > .01 
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Equation 

No. 


6 . 


7. 


y. 


10 . 


i l . 
12. 


13.. 


Equation 


Test Option Flag C 
If Flag C = 0, Pump 

Go to Equation 13. 

If Flag C = 1, Collector 
Continue. 


N 1 * /q 
(npsh) 3/4 


fc _ . f (s Q V ) 

NPSH f 2 (S ’ N x } 

(See Curve 67-135 as example) 


AH - 1< 2 (NPSH) f 


AH. 


NC 


N„ 


*2 


f ( _ a _ jl 

3 N., * * N, 


N 1 "“ “ “l 1 

(See Curve 67-132) 
ah N c = k 3 n / 2 f 3 


AH - AH nc - ;VH C 


P = P + WJL 

*<1 SH 144 


T = T 
e J S J 

EXIT 


Main Pump Equations 
Using first past time value of AH^. 
^ ? AH* 


Q = K (AH*) 1/2 
r 

AQ/At =CQ (t _ 1) - Q ((; _ 2) yAt 


i'f. 


AQ 

5 


*•' , r l JL 

i At’ Q 


-) 


nom 

■f.-iyt’** i.iu’.’p 67 -3 31 as 


example) 


Boundary Conditions 
or Limits 


Min 



< Max 


These will be 
negative if collector 


p > p 
d - V 
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Equation 


Boundary Conditions 
or Limits 


1 w Re 

Q = Q [1 + K x — — ] + Q j 

(NPSH) 3 ^ 4 

= f (s d 

NPSH 2 ' * N 1 ’ N x 

(Curve 67-135) 

AH C = K 2 (NPSH) f 2 

^NC _ . £ 

2 3 'N * N 

N 4 J "l 1 

(See Curve 67-133) 

4H N C - K 3 H 2 l 


If Q <0, error 


0 < S < Max 


0 < < Max 

""i - 

t} 2 * 

0 < — — < Max 
_N i - 


Test Option Flag B 
If Flag B = 0 

Go to Equation 21. 

If Flag B - + 

Continue. 

AN.. /At » (N - N. 

(t) X (t-1) 


= f ( _ jj -2—) 

„ 2 *4 'At ’ 1’ N/ 

1 

(Optional Curve) 

AH = K. N 3 f. 
tr 4 14 


( <n V _aL 
4H W 5 V N 1 ’ W 

(See Curve 67-125 as example) 


-1 < k - < 4-1 

- 54H kc- 


4I W - 4 V (1 - K 5 
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Equation 

No. 

23. 

24. 

25. 

Zb , 

27. 

28. 

29. 

30. 


Equation 

4H - AH^l - AH, - AH tr 


Boundary Conditions 
or Limits 


P = P + M-EL 
a r SM 144 


P , > P 
d — 


V 


Before w a Convergence: Bypass to waterhammer pump. 

After w a Convergence: Continue. 


1 N * 

x — f re 

p NPSH 6 v ’ N ’ N ' 
(Curve 67-136) 

Ap c = K 6 p (NPSH) f & 


P N -, 


, V 

f 7 ( N 1 ’ 


N„ 


(See Curve 67-134) 
p 1 = K ? p N ; 2 f ? 


Test Option Flag B 
If Flag B = 0 

Go to Equation 29. 
If Flag B = + 

Continue. 


Ay. N.. 

-7i * f 8 <-r> “i- % 

pN. 1 


(Optional Curve) 

Ay = K. p N 2 f 
tr x 1 


8 


p = y - Ay c - Ay fcr 


Test Option En Flag 
If H 2 Flag = 0, Incompressible 
Go to Equation 32. 

If H 2 Flag = 1, Compressible 
Continue . 
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Equation 

No. Equation 

31. COMPRESSIBLE: Calculate T ,.Ene, p 

Method is the same as outlined under 
Inducer Subroutine, except 

y (3()) N 1 

AE n u = T78”w~" 

EXIT 


32 „ 


INCOMPRESSIBLE: Calculate T , p 

e e 


T + — 
s C T 


N 1 

1 778 w 


= f 


(T ) 
v e 


AH . 
77 8 J 


EXIT 


Boundary Conditions 
or Limits 


T > T 
e — o 
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FEED SYSTEM SUBROUTINES 
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I. GENERAL ASSUMPTIONS 


A. Constant density in all lines. 

B. Constant area in each section. 

C. Pressure waves travel at sonic velocity (i.e. , the flow velocity 
is negligible compared to sonic velocity) . 

D. Constant elevation and length. 

E. Conduit walls are elastic laterally but not axially. 

(No Poisson's ratio effects.) 

F. Fixed computing internal, with wave propagation times for each 
section being an integral number of this interval. 

G. Static pressure equal to total pressure (for simplified presen- 
tation of equations only, program includes effect) . 

H . No friction in line, instead it is lumped at the joints. 

II. BASIC EQUATIONS 

Using the above assumptions, the basic equations for a simple elastic- 
conduit are 


A 

8 1 

f T~r T / T 

/ / 7 7 7 7 7 / ""7"7*1 

4 V* 

— 1 1. 

1 1 

CZ3ZZ=Z_J7^ 

-Z... / / / / /...■/. V / Q-J 


- * — *> X 

+ — Z 


E'P’vtion of Miction: 


aH _ _ i_ av 

3X ~ g 9t 
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Equation of Continuity: 


Where : 


w 3V 
g 8X ’ 


(3) V 


,, , Bk.d \ 

p (1 + "ITth > 


By methods suggested by L. Bergeron, "Waterhammer in Hydraulics and Waves 
Surges in Electricity" translated by AsME, 1961, Wiley, New York, these 
equations can be solved in a finite difference form 


(AN - — a v . it = — , (5) s - — 

w At At w’ V ’ ^ ' gA 

w 


for right traveling waves 


t t + ^A,t* W A t* 

(6) + S -2-i^CC = • + S > • 

P P p p 


for left traveling waves 


P w P w 

(7) -Alt . s = c- - -i£S. - s -2*£i 

P P P P 


where 


t = current time 


t* = (t - At) « t - 


Equations (6) or (7) can be solved with a known boundary condition or another 
equation relating Px, t and wx,t» 

A computing interval, (At) m £ n , is selected which when multiplied by an 
interger (1^) will equal the At of each section, 

L i 

At- = — — , I- (At) . “ At, 

IV’ 1 mxn 1 

ws 


it - — , i «t) - it 2 

wd 


At 3 " V * I 3 (At) min At 3 
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III. 


DESC RIPTION OF COMPONENTS 
A. TANK AND INLET LINE, SR NO. 1 

ft 


-4 


• ** 




>v, 




L 


3 


< 1 



1 Pa 

f 

> 

- ~ VlO 


Pb 


1. Boundary Conditions 

U) (i) PT = PA or (ii) P 


•a 2 

P - ™ R 
A T p 


(2) Known pant time conditions of P , wB 

B 

Initial Conditions 


(1) w A 0 

(2) w B * 0 


P = P = P 
A B T 


(i) ? A “ P.J, calculate w^, C ' 


3. Equations 

(1) t* = t 


w 


12) C, 


_T. _ g wA, t 
P ° P 


<3) c ' - -£«i£ - s iaiai 

A P P 


Solving for wA, t from eq (2) 

“» » A .t - f $ ■ '•„-) 


Note: the. si xn of (~ C, ) determined the direction 

i • A 

;?t the. flow, Plow out- of tank is assume '* 

:)OS » t ;i ve . 
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(ii) P A - 

v 

P^ - — R (Friction loss at tank exit) 

4. Equations 

(1) 

* L 

t - t - ~ 

(2) 

P 

P A,t P T p R 

(3) 

-Ail - s * c - 

P p A 

(4) 

P W 

c ~ - s 

A p p 


Solving equations (2) and (3) for w , gives the solution to the quadratic 

A, t 

(5) { - 1 + Z 1 + 71 ( r - c a']i 

Because friction represents a loss the positive root is selected. 

Also, the direction of the flow is determined by 

(6) w A>t - sign{£ - c A -)} |i A(t | 


as was shown previously. 


Then, 


w 


(7) P A,t ' » (C A _ + S 




B. INTERNAL POINT WITH OR WITHOUT FRICTION, SR NO. 2 


4*<3 "" — A* 
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1. 


2 . 


Known Conditions 

a. Without Friction 

(1) w B - w c 

(2) P = P (velocity head neglected) 

B c 

(3) Known past time values at A and D 

b. With Friction 
(D w B * w c 

w 2 

(2 > P C ■ P B - P R 
(3) Known past time values at A and D 
Equations 


a. Without Friction 


(1) 

V 

4 . AB 

" V . 
wA 

(2) 

t* 

_ Lcd 

B 

VwD 


P t 

w 

(3) 

B, t 

- -t- s — 5 jJ 


P 

A p 


C + = + s 

A p A p 


P w P 

(A) Sul _ s -Jjt. = c " = 

^ } p D p D p 


c Vt*D 

0 p 


(5) P " P 
K J B,t c,t 


■i ' s 1 (2), (4), and (5) for gives 

o > c 


W B,t = P 


: A + - V 

SA + SD 

•* 


( 6 ) 



Appendix B 


b. With Friction 

(1) - (4) are the same. 

(5) p c- p b- ! t r 

Solving eq's (3), (4), and (5) for w R gives 


(6> K.tl-'ai <S A + V{ - 1 + J 1 + l^> r |(C A + - C d"} 

the sign or direction of the flow is determined by the difference between 
the characteristics. 


(7) SIGN = SIGN 


therefore 


K + - c d'>} 


(8) * B>( . = SIGN | w B>t 


then 


w r 


(9) P B - P (C/ - S A -|*£) 


W T 


do) p c - p (c„- + s D -£^; 


C. BRANCH POINT, SR NO. 3 
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(a) Equations 


U) w B = w E + w F 

(2) P = P^ = P_ (change in Kinetic energy neglected) 
a b E 


(3) t A * = t - v 


AB 

wA 


(4) t * «■ t — — 

V ' D V V _ 
wD 


(5) t * ■= t - 

wc 




P A * 


B > . *• + s >5 = c + = + s —^ JL 


w * 
A.tA 


A P 


A p 


P W p * w * 

(7 ) 4^- S n JSat „ c - = s. 


D p D 


D p 


(8) 


P w p * 

_I*t . s Jjt . , c - 


w „ * 

c ,tc _ s c ,tc 


p c p c p c p 

Solving (1), (6), (7) and (8) simultaneously 


< 9 > V p 


c + C c 
_L. + _JL + _£ 

SA SD Sc 
SA SD Sc 


then, 


C A " P B t /p 

(10) *B,t = - SA/p' (from eq (6)) 
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K t ■ - n 

(11) w E>(; = ( — J C D ) (from eq (7)) 

(12) w F>t - w B>t = w Ejt (from eq ( 7 )) 

Note: The equations as presented are in the form used in the 

computer. This was determined to be convenient form. 
It is recognized that some of the equations could be 
simplified by algebraic manipulations. 

Also Note: The units of pressure (P) in these equations are 

lb/ft^. In the program the units are lb/in. 2. 

therefore 

P (eq'ns) ” P (program) X ~2 
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D* CONTROL VALVE, SR NO. 8 



a. Equations 


(1) W, 

(2) w. 


B 


w 


B,t 




62.4 


(3) K 


w 


K *.(t), (K is a known function of time) 
w w ■ 


(4) R 


62.4 

K 2 
w 


(5) 

( 6 ) 

(7) 

( 8 ) 




P 


c,t 


) 


V = 


V = 


- ( 12 ) 


t 


t 


AB 

^WA 

l cd 


WD 


These equations are the same as those listed in, 
SR it 2, (b)(ii). Equations with friction, 
Equations (6), (7), (8), (9), and (10). 
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WATER-HAMMER PUMP, SR NO. 11 



a. Known Conditions 

(1) ■» N(t) (known pump speed or speeds) 

(2) T = K (constant inlet 


(3) w. 


B 


= w - 0 
c 


b. Equations 


(1) w„ = = w , P D „ = P_, P„ „ * P. 


( 2 ) 


B c B,t S c,t 

L, 


t * = t - 

A V 


AB 

wA 


(3) t D * - t - v 


CD 

wD 


< 4 > p s W P P A t A W A tA* 

— + s — = c + = A>tA + s A,tA ■ 

p A o A p A p 


(5) Call NASA-TP 


P d ' P S + “l- V T o > 


The head rise in the pump can be calculated 
as a function of wp t with known value of 
speed, temperature, etc. 
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(6) P 


C Zp - r ~ = * D » tD * _ q ISjtD* 
S D p “ C D p b D p 




Equations (4), (5) , and (6) t can be solved 
simultaneously for P g , P^, w p by internal 
iteration as follows. 

i) Assume w_> 


ii) From equation (4) calculate Pg 

iii) Call NASA-TP Subroutine, and calculate P^ 

iv) Using calculated value of P^, calculate w p 
from Equation (6) 


• P 

■n n 


(rr ~ Cl) 


“D P 

v) Compare calculated value of w p with assumed value until 
W P W P (assumed) e 

vi) Call Iteration 2, the convergence subroutine and 

iterate until the allowable error (&) is less than 
or equal to 

e = Max (| e| , x w p ) 
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F. OUTLET SR NO. 12 


+VL 




L-a© 

“T 

5 


S' 


?e. 


i 

ft 


I 

& 


1. Boundary Conditions 


(1) P « P (known exit or atmospheric pressure) 

6 3 . 

(2) Known hydraulic resistance, R (friction or orifice) 
2. Equations 

a. Without Friction 


(1) t* - t - 


AB 

wA 


V 


(2) ISlJL + s, it . c + - + S, ^ 

p A p A P A p 


(3) P B t - MAX (P v> P & ) (P v = vapor pressure) 

w »B, t - 1; (0 a* - 


b. With Friction 


(1) and (2) are the same 


(3) P, - MAX(P , P„) + 




B,t 


(4) jw. 


B , t 


P S A 

2R 




4R 


c a + - -J 
A p 
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(5) 

SIGN » 

SIGN 1 

i 

k * 

-ij 

P ^ 

(6) 

”B,t S 

SIGN 

K, t 

1 

(7) 

P B,t * 

0(C A 

- S H 

it, 

p 


NOTE: 

All points are tested in the program 
for vapor pressure 


TEST (P. „ - P ) 
i,t v' 

If the pressure is less than P^, it is set equal to P 


h.t = «“ V 

and all equations using it are re-evaluated. 




APPENDIX C 

LISTING OF PROGRAM DECK 
AND 

TEST CASE DATA 
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I M P L I C I T W h A L * B < A - H , f ) - / ) 000 1 ( ) 

CflMMON/ARFA 1 /DA TAG ( 1 0 0 ) ,C JRVFS ( 1 000 ) , DATA ( 1 ?00 0 ) 0 00 >0 

COMMON/ ARE A? /T I T L Hi ( 1 6 ) ,LPRT ( 200 ) ♦PPTT ( 100 ) 00030 

1 LD ( 1 00 ) , N1 ( 1 00 > *'M1 ( 1 00 ) * MPR ( 100 ) * NMF ( 100 ) 0 0 040 

COMMON/ ARHAWH/RHO ♦RH01 44 * L ST ART ♦ L JO I NT , NS ♦ LOC 000S0 


*.r 

vC 

c 

G 

C 

C 


CCCCCCOOCCC 
C O C O r-i ,-1 r- 

occcco. coocc 
c c c g o c c o c c e 


VC 


o c 

< I 

t— Q 
<£ 4- 

m d * 

a; + + 

\a m 

N 

♦ in 
•o c 

a < ~ 

~ f c o 

< < c 
<- a ..'nacFH 

a o m 5 ii < n 

~ f in co m n 4 m 

a a ijj < t h h a <t 

.(A 1 I- h 3 ^ 'H i/? 

< C ~ w C 

-J J J D -J < < I- J H 

<<<n<4<a<cz 
UUU~UOD.J<JC 'Ii 

in in o in in o it in o 

a, <t in m n x x a o 
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oooooooooooooooo o 

NflDO'OH(Mr 0 4 inONflDO'O^RlfO 
^ HHWwtvMcywwMMMnwfon 
o oooooooooooooooo 
c oooooooooooooooo 


0000000.0 00 o'o o o o o 

«V)ONtt&O'^CMm«tn.iO.N<D 0 t 
ro ro ro 10 n ro 4 4 < * 4 <t 

oooooooooooooooo 
oooooooooooooooo 


OOOOOOOOOOOOQP 

OrtW« 4 K)' 0 K®^O»-flJI 0 

lOjninininininintfiino^oo 

OQQOOOOOOQOOOO 

OOOOOOQOOOOOOO 



MPfiM)=0 00640 

NWPiI1=0 00650 

102 LPRT* I*100)=0 00660 



DO 103 1=1*1000 00670 

103 CURVES* I >-0* 006H0 

N = N(JC 0 MS #100 +1 10 0 0 690 


c 

c 

c 

c 

o 

C 

c 

c 

C 

c 

o 

c 

o 

c 

c 

c 

c 

o 

c 

o 

c 

c 

© 

c 

c 

p 

c 

© 

o 

c 

© 

p 

c 

c 

c 

c 

p 

c 

c 

© 

© 

c 

c 

© 

C 

c 

c 

PM 

(V 

x 

vt 

IT 

x 

a 

00 

O' 

c 

PM 

CVi 

x 

<t 

X 

X 

N 

X 

O' 

© 

a 

<v 

ro 

■St 

10 

X 

N 

X 

0 

© 

•Pi 

CV 

ro 

< 

X 

X 

a 

X 

o 

© 

PM 

CV 

r r - 

<3 

if 

IS 

a 

N 

a 

N 

IS 

N 

N 

IS 

N 

00 

cc 

00 

oc 

oe 

oe 

CC 

OC 

cc 

GC 

O' 

O' 

O' 

O' 

O' 

o 

0 

o 

o 

O' 

c 

C 

o 

c 

© 

c 

c 

c 

c 

c 

— * 

— 

PM 

p~ 

.PM 

PM 

o 

© 

c 

o 

O 

o 

C 

O 

C 

O 

o 

o 

o 

o 

o 

c 

© 

O 

o 

o 

C 

o 

© 

o 

© 

Q 

© o 

o 

p 

..**• 

r~> 

PM 

rH 

pH 


"H 

«-* 


•— 

>— 


— 1 

PM 

r-* 

pp 

© 

c 

c 

6 

O 

o 

c 

C 

C 

O 

o 

c 

o 

o 

c 

o 

C 

o 

c 

a 

© 

© 

© 

c 

o 

o 

c 

p 

p 

c 

o 

o 

© 

© 

c 

:C 

c 

e 

c 

c 

c 

c 

c 

Q 

c 

o 


£V 

+ 

y 

X- 

UJ 

> 

a 


> 

a 


z 

c 

a 

LL 

C 

< 

lii 

a 

w 

< 


© 

CV 


_i 

a 


t r. 
a 
> 
a 


y ~ 
— 

x + 

nr. ^ 

CV ~ 






*> 






• 

PH 










c 


© 




CV 

• 

eg 



MM 





li 

X. 

• 

c 








c 






#> 

HP 








© 


X 


O' 




» 

PM 

» 



a 





a 

u 

y 

X 








CVJ 

pH 





CV 

a 










PH 


pH 




X 

.MM 

X 



PH 




G 

a 

> 


eg 








PM 






MM 

a 

a 









» 


•> 



PM 

pH 

HP- 

CV 



M 




X 

a 

a 

- 

• 








M 

| 




h* 

•> 

a 

cr 







.#• 


n 


X 



pH 

CV 

a 

CV 



CVi 


* 


CV 

*> 

z 

H 

c 








pH 





a 

• 

_j 

< 







m 


X 


O' 



+ 

* 

_j 

• 





-H- 


r 

X 

u 

PM 

<t 








) 

PM 




< 



H 


PH 

PH 

PH 

PH 

PH 

N 

in 

«H 

X 

pH 

X 


PM 

X 

a 

X 

PM 

PH 

• 


•a 


X 

X 

* 

P 

cv 








PH, 

H 

MM 



K 

Q 

PH 

to 


pH 

CVJ 

ro 

<t 

in 

pH 

> 

» 

t— 

* 2 


HP 

PH 

MM 

CV 


X 

X 

PH 

* 

c 

X 

cv 


X 

r 








a 

tr 




X 

a 

o 

z 


CV 

CVJ 

CVJ 

CV] 

CVJ 


a 

m 

a 

X 

o 


z 

eg 

a 

CV 

CV 

CV 

W 

X 


CV 

CV 

1! 

X 

V 

o 








a 

< 

*■ 



z 

o 

» 



HP 

HP 

w 

HP 

HP 

in 

u 

X 

a 

a 

u 


u 

p** 

HP 

PM 

HP 

w 

a 

< 

X 

► 

p 

c 

+ 

a 

X 



y 

z 




< 

a 

*H 




u 

h 

If 


u 

u 

u 

u 

u 

N 

c 

vi 

Q 

»H 

o 

X 

It 

PH 

PM 

PM 

u 

u 

a 

X 

a 

pH 

X 

z 

y 

a 

cv 



PH 

*• 




H 

X 



x 

pH 

UJ 

X 



2 

z 

z 

2 

2 

pH 

z 


z 


z 

• 

PH 

O' 

X 

c 

2 

2 

H 

PM 

> 


X 

a 

HP 

*M 

PH 



a 

pM 

• 



X 

z 

It 


,pH 

CVJ 

a 

• 

PH 


II 

II 

ii 

II 

11 

p» 

M 

**+ 

it 

PM 

II 

.pH 

X 

X 

*M 

O' 

II 

II 

a 


a 

II 

cv 

•• 

X 

a 

X 



X 

1) 

c 

p 

ro 

z 

V 


PH. 

d 

PH 

© 

X 

pH 

© 

PH 

PH 


PH 

PH 

in 

PM. 

X 

-s 

X 

PH 

II 

PM 

PM 

» 

pH 

p*> 

PM 
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